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SUMMARY 

Results  of  transonic  w ind  tunnel  tests  on  two-dimensional  aerofoils  are  analysed  with 
respect  to  the  effects  of  wind  tunnel  wall  interference.  The  tests  were  conducted  on  two 
geometrically  similar  models  of  each  of  two  aerofoil  sections  the  NACA  0012  and  the 
RCK-I  sections  and  covered  a range  of  Mach  numbers  from  0-5  to  0-82  with  model 
chord  to  tunnel  height  ratios  of  0 • 125  and  0-25.  Results  from  measurements  on  all  models 
in  both  solid-  and  slotted-wall  test  sections  are  corrected  for  wind  tunnel  wall  interference 
effects  by  the  application  of  classical  linearized  theory.  For  the  solid  wall  results,  these 
corrections  appear  to  produce  data  which  are  very  dose  to  being  free  of  the  effects  of 
interference.  In  the  case  of  slotted  walls  however,  linear  theory  is  found  to  significantly 
underestimate  the  magnitude  of  blockage  interference.  The  introduction  of  an  empirical 
blockage  correction,  when  combined  with  the  linear  theory  representation  of  lift- 
interference,  enables  the  cross-flow  characteristics  of  the  slotted  walls  to  be  determined, 
and  leads  to  a successful  correction  scheme  for  the  slotted  wall  results.  Extensive  com- 
parisons with  data  from  other  sources,  both  experimental  and  theoretical,  provide  further 
verification  that  the  corrected  results  are  indeed  closely  interference  free. 
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NOTATION 

Area  of  aerofoil  section  in  v ; plane 
I fleet  is  e area  of  aerofoil  section  in  v ’ plane 
Width  of  sloi  (see  Fig.  3) 

Breadth  of  tunnel 

Wake  drag  coefficient  (drag  force  from  wake  traverse)  } pU*S 

Drag  coefficient  (drag  force)  j pi^S 

l ift  eoelticient  (lift  force)  J pi  -S 

Normal  force  eoelticient  (normal  force)  \pU*S 

Pitching  moment  eoelticient  (pitching  moment  about  <■  4)  J pU'Sc 

Pressure  eoelticient 

Aerofoil  chord 

Constants 

Mean  periodic  spacing  of  slots  (see  big.  3) 

Universal  empirical  blockage  factor 

Slotted  tunnel  geometry  factor 

Ratio  of  corrected  to  uncorrected  kinetic  pressures 

Free  stream  total  pressure 

Height  of  tunnel 

Integrals  (see  Section  3.2) 

Constants 

Aerodynamic  loading  ( /»/  pu)  \pU2 

Lift  per  unit  span  of  aerofoil 
Thickness  of  slats 
Free  stream  Mach  number 

Strength  of  doublet  representing  aerofoil  displacement 

Strength  of  source  representing  aerofoil  wake 

Effective  number  of  full  width  slots 

Number  of  measured  values  of  Mach  number 

Wall  porosity  parameter 

Local  static  pressure 

Free  stream  static  pressure 

Dummy  integration  variable 

Reynolds  number  based  on  <■ 
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\k  mtt  area  (/>  ■ c) 

I rce  stream  temperature 
Aerofoil  maximum  thickness 
f ree  stream  velocitv 

St  ream  vsi  sc  velocitv  induced  b>  tunnel  walls 
Vertical  velocitv  induced  b>  tunnel  walls 
Streamwise  co-ordinate  of  tunnel,  origin  at  mid  chord 
Chordwise  co-ordinate,  origin  at  leading  edge 
Spanwise  ordinate,  origin  at  mid  span 

Aerofoil  thickness  co-ordinate,  also  vertical  co-ordinate  of  tunnel 
Angle  of  incidence 

I’randtl  Cilauert  compressihihtv  parameter  (1  \l  ■)' 

Circulation 

Ratio  of  specific  heats,  taken  to  he  I 4 

Prefix  denoting  correction  increment 

Root  mean  square  deviation  from  mean  curves 

I ill  interference  parameter  associated  with  stream  directum 

I aft  interference  parameter  associated  with  streamline  curvature 

Blockage  factor 

transformed  co-ordinate  (see  Section  3.2) 

Chordwise  parameter  tsec  Section  3 3) 

Doublet  strength 

Transformed  co-ordinate  (see  Section  3.2) 
t ree  stream  densitv 
Perturbation  velocitv  potential 

Perturbation  velocitv  potential  due  to  aerofoil  in  unconstrained  flow 
Perturbation  velocitv  potential  induced  bv  tunnel-wall  constraints 
Ratio  of  ventilated-wall  to  closed-wall  values  of  blockage  factors 


Subscripts 

H Denotes  total-blockage  values 

c IVnotes  closed -(solid-)  wall  values 

f Denotes  free-air  (i  c corrected)  values 

I Denotes  aerofoil  lower-surface  values 

« Ifenotes  solid-blockage  values 

u Denotes  aerofoil  upper-surface  values 

»•  Denotes  ventilaied-(slottcd-)  wall  values 

m IVnotcs  wake  blockage  values 


I.  INTROIHC 'HON 


I he  derivation  of  two-dimensional  aerofoil  characteristics  from  wind  tunnel  tests,  especially 
at  transonic  speeds,  is  suhiect  to  several  sources  of  experimental  error,  one  of  the  most  significant 
of  which  is  wind  tunnel  wall  interference  When  an  aerofoil  is  placed  in  a wind  tunnel,  forces 
and  pressures  measured  on  it  differ  from  those  which  would  lie  obtained  if  the  model  were 
tested  under  similar  conditions  in  free  an  because  the  wind  tunnel  walls  interfere  with  the  flow 
It  is  usual  to  make  "corrections"  to  the  measured  quantities,  and  to  Mach  number  and  angle 
of  incidence,  so  that  the  "corrected"  results  might  correspond  to  those  which  would  bo  obtained 
if  the  model  were  tested  in  free  an  at  the  "corrected"  Mach  number  and  angle  of  incidence. 
Such  corrections  mas  be  reduced  by  using  models  with  chords  very  much  less  than  the  tunnel 
test  section  height  However,  in  an  otfori  to  make  maximum  possible  use  of  the  available  tunnel 
Reynolds  number,  models  are  generally  designed  with  large  chord  to  tunnel  height  ratios,  giving 
rise  to  significant  errors  due  to  wind  tunnel  wall  interference 

I he  problems  in  obtaining  accurate  two-dimensional  transonic  aerofoil  data  from  wind 
tunnel  tests,  and  m attempting  to  eliminate  uncertainties  in  transonic  two-dimensional  aerofoil 
theory  were  recognised  by  the  Commonwealth  Advisory  Aeronautical  Research  Council.*  1 hat 
body  endorsed1  a proposal  by  the  Co-ordinators  in  Aerodynamics-  to  "test  a small  number 
of  suitable  aerofoils  in  as  many  as  possible  of  the  transonic  wind  tunnels,  throughout  the 
Commonwealth,  in  which  two-dimensional  research  is  normally  done".  \ proposal  for  such  a 
research  program,  put  forward  bv  the  United  kingdom  CAARC  Co-ordinator  in  Aerody  nannes,3 
stated  that  the  purpose  of  such  a program  would  be  "to  add  to  our  mutual  knowledge  and 
understanding  of  experimental  techniques  in  a subject  which  is  at  once  of  great  practical 
importance  and  at  the  same  time  afflicted  with  serious  uncertainties,  both  theoretical  and 
experimental". 

I wo  basic  aerofoil  shapes  were  selected  for  the  co-operative  program.  I he  first,  the  NACA- 
(X>I2  section  was  selected  to  provide  a "standard"  symmetrical  section  I he  second,  an  aerofoil 
designed  bv  the  numerical  hodograph  plane  method  of  Bauer,  Ciarabedian  and  Korn,4  is  an 
example  of  a modern  "supercritical"  section,  with  a reasonable  extent  of  shock  free  supersonic 
flow  on  the  upper  surface  and  some  degree  of  rear  loading  I his  aerofoil  (hereafter  referred  to 
as  Blik-I)  had  the  added  advantage  of  having  already  been  extensively  tested  at  high  Reynolds 
number  in  the  NAI  wind  tunnel  in  Canada  A*  Since  both  these  sections  already  formed  part  of 
a theoretical  and  experimental  investigation  of  transonic  aerofoils  at  ARl  . with  models  of 
each  either  completed  or  being  manufactured,  it  was  decided  to  extend  the  experimental  investi- 
gation to  include  the  test  program  proposed  m Reference  .V  Results  from  tests  in  a slotted 
wall  tunnel  on  two  models  of  differing  chord  lengths  of  aerofoil  BCik-l7!tM  and  of  three  models 
of  the  NACA-001  ’ section1"  have  been  reported  previously  These  tests  were  later  repeated  with 
the  slotted  walls  replaced  with  solid  ones  11 

The  present  report  concentrates  on  the  analysis  of  the  above  test  results  with  respect  to 
wind  tunnel  wall  interference.  Tor  a model  of  fairly  small  chord  to  tunnel  height  ratio,  the 
corrections  can  be  calculated  theoretically  by  linearized  potential  flow  theory,"  in  which  the 
changes  to  the  stream  direction  and  the  streamline  curvature  are  referred  to  as  lift  interference, 
and  changes  in  stream  velocity  and  longitudinal  gradient  are  referred  to  as  blockage  interfer- 
ence The  lift  interference  arises  from  changes  induced  by  the  tunnel  walls  in  the  circulation  or 
vorticity  around  the  model  and  is  assumed  to  be  independent  of  the  blockage  effect,  which  is 
associated  with  changes  in  the  velocity  potential  of  a doublet  and  source  representing  the  volume 
occupied  by  the  aerofoil  and  its  wake 

A major  problem  to  tv  faced  when  applying  such  linear  theory  to  a practical  wind  tunnel 

* Australia,  Canada,  India.  New  Zealand  and  the  United  kingdom  arc  active  members  of 

CAARC. 
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situation  is  that  of  correctly  representing  the  boundary  conditions  at  the  tunnel  walls.  For 
the  case  of  solid  walls,  the  boundary  conditions  are  straightforward  and  mathematically  amen- 
able. However,  in  transonic  flows,  porous  or  slotted  (as  is  the  case  for  the  ARL  transonic  wind 
tunnel)  walls  are  often  used  Such  walls  are  usually  arranged  primarily  for  the  generation  of 
low  supersonic  Mach  numbers  using  diffuser  suction,  or  to  minimize  shock  reflections  at  high 
subsonic  and  low  supersonic  Mach  numbers,  with  little  attempt  being  made  to  utilize  the  avail- 
able reduction  in  inteference  made  possible  by  the  non-solid  walls.  In  this  case  the  wall  boundary 
conditions  are  complicated  and  it  is  usual  to  change  the  mixed  boundary  conditions  into  an 
equivalent  homogeneous  boundary  condition  which  produces  the  same  solution  at  the  centre 
of  the  test  section,  but  the  correct  representation  remains  open  to  question. 13-13  As  data  are 
available  from  tests  in  solid  as  well  as  in  slotted  walls,  it  should  be  possible  to  examine  the 
validity  of  the  assumption  of  the  independence  of  lift  and  blockage  interference,  and  hence  the 
applicability  of  linear  theory,  separately  from  the  physical  and  mathematical  representation  of 
the  non-solid  wall  boundary  condition,  since  the  solid-wall  boundary  condition  is  straight- 
forward. 

The  statement  of  the  results  of  linear  interference  theory  given  in  the  standard  references 
has  until  recently  included  errors,  originating  in  the  early  analyses  and  faithfully  reproduced 
by  later  authors.  These  errors  were  pointed  out  by  Catherall.11  It  was  therefore  considered 
worthwhile  to  include  a rev  iew  of  the  major  results  of  classical  linear  theory,  and  this  is  included 
in  Section  3. 

2.  EXPERIMENTAL  DETAILS 
2. 1 Section  Definition 

Aerofoil  section  BGK-1  was  one  of  the  earliest  designed  by  Garabedian  and  his  colleagues 
using  their  numerical  hodograph  method.4  The  section  was  designed  to  produce  shock-free 
flow  on  the  upper  surface  at  A 1 = 0-75  and  <*  = 0°,  giving  a Cl  of  0-63  ( i n viscid ) (Fig.  1).  The 
ordinates  of  the  aerofoil  were  derived  from  those  published  in  Reference  5,  with  the  exception 
of  those  close  to  the  trailing  edge,  where  the  presence  of  a cusp  required  some  modification  to 
simplify  model  manufacture.  These  modifications  involved  thickening  the  section  in  the  region 
0 95  < atl/c  ^ I 00  and  are  noted  in  Reference  9. 

The  NACA-0012  is  a symmetrical  section,  12%  thick,  with  ordinates  given  by: 

* = 0 6 jo-29690^ -0  12600  *L  - 0 35160  (^jJ  + 0-28430  -0- 10150 

with  a leading  edge  radius  of  0 0158c.  The  basic  section  has  a blunt  base  of  0 00252c.  However, 
to  avoid  confusing  comparisons  between  models  arising  from  base  pressure  variations  with 
Reynolds  number,  the  upper  and  lower  surfaces  were  linearly  extended  to  a sharp  (01  mm 
thick)  trailing  edge.  Due  to  this  extension  the  actual  physical  chord  length  of  the  models  was 
slightly  greater  than  the  nominal  chord.  The  calculated  force  coefficients  are  based  on  the  nominal 
chord  (c). 


2.2  Model  Dimensions 

The  testing  of  two-dimensional  aerofoils  in  transonic  wind  tunnels  presents  conflicting 
requirements  on  model  chord;  a small  chord,  while  minimizing  tunnel  interference  effects, 
introduces  difficulty  in  producing  a model  of  sufficient  strength,  and  limits  the  number  of 
pressure  orifices  which  can  be  included.  A larger  chord  also  allows  greater  use  to  be  made  of 
available  tunnel  Reynolds  number.  Model  aspect  ratios  should  also  be  large  enough  to  produce 
sensibly  two-dimensional  flow  conditions,  at  least  over  the  central  measuring  section.  Since 
one  of  the  aims  of  these  tests  was  to  investigate  the  effects  of  tunnel  interference,  models  with 
larger  values  of  model  chord  to  tunnel  height  ratios  than  would  normally  be  tested  were  con- 
sidered desirable,  both  to  increase  the  magnitude  of  any  interference  making  it  easier  to  measure, 
and  to  investigate  the  limits  of  applicability  of  linearized  theory.  On  the  basis  of  these  con- 
siderations, a maximum  model  chord  of  203  -2  mm  was  selected.  The  cjh  ratio  of  0-25  is  within 
the  suggested  limits  of  applicability  of  linearized  theory  (cjh  < 0-3), 3 and  the  aspect  ratio  (2-6) 
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is  sufficiently  large  to  avoid  major  effects  of  three-dimensionality.  Models  of  both  sections 
with  this  maximum  chord  were  therefore  tested,  as  were  models  with  101  6 mm  chords 
(giving  < h 0 125,  aspect  ratio  5 • 25)  to  provide  comparative  results. 

2.3  \N  ind  Tunnel 

All  tests  were  carried  out  in  the  ARL  variable  pressure  transonic  wind  tunnel.  The  test 
section  fitted  for  the  slotted  wall  tests  had  solid  side  walls  and  longitudinally  slotted  top  and 
bottom  walls  (open  area  ratio  16  5"„  at  the  model  location),  with  dimensions  at  the  model 
location  of  h 533  mm.  h XI 3 mm  (Fig.  2),  surrounded  by  a plenum  chamber  2540  mm  in 
diameter.  These  walls  are  those  normally  used  for  two-dimensional  testing  in  this  tunnel.  The 
open  area  ratio,  which  was  selected  primarily  for  the  generation  of  low  supersonic  Mach  numbers 
through  diffuser  suction,  is  known  to  be  larger  than  that  giving  interference  free  subsonic  test 
conditions.  For  the  slotted  wall  tests,  Mach  number  was  derived  from  measurements  of  the 
pressure  in  (he  plenum  chamber  surrounding  the  test  section,  and  in  the  entry  to  the  contraction, 
assuming  these  to  be  the  static  and  total  pressures  of  the  test  section  flow  respectively.  For  the 
solid  wall  tests,  test  section  total  pressure  was  once  again  derived  from  the  static  pressure  in  the 
contraction  entry,  while  test  section  static  pressure  was  measured  at  a side  wall  pressure  hole.14 

2.4  Transition  Fixing 

The  effects  of  transition  fixing  on  the  BGK-I  section  have  been  studied  previously.8  That 
investigation  showed  the  need  to  fix  transition  if  gross,  shock  induced,  laminar  separations  were 
to  be  avoided  at  Mach  numbers  and  angles  of  incidence  above  the  design  conditions  (A/  0-75, 

» 0 ).  However,  the  use  of  transition  fixing  destroyed  agreement  between  theory  and  experi- 

ment at  the  design  condition.  In  spite  of  this  problem,  it  was  decided  to  conduct  the  slotted  wall 
tests  with  transition  fixed  on  both  surfaces  to  avoid  both  shock  induced  laminar  separations, 
and  the  misleading  variations  of  drag  and  pitching  moment  (and  to  a lesser  extent  lift)  caused 
by  the  movement  of  the  chordwise  position  of  transition  with  Mach  number  and  angle  of 
incidence.  Also,  the  transition  free  behaviour  of  the  section  had  already  been  the  subject  of  a 
detailed  investigation.7 

Transition  was  therefore  fixed  on  both  surfaces  of  the  BGK-I  models  for  the  slotted  wall 
tests.  For  consistancy,  a similar  transition  fixing  was  used  for  the  slotted  wall  tests  on  the 
NAC'A-0012  section,10 and  these  conditions  were  maintained  throughout  the  tests  in  solid  walls.14 
Surface  oil  flow  observations  indicated  that  at  the  maximum  lest  Reynolds  number  for  each 
model,  boundary-layer  transition  was  effectively  produced. 


.3.  SI  BSON1C  LINEAR  INTERFERENCE  THEORY 
3.1  Introduction 

In  subsonic  linear  theory,  interference  is  assumed  to  be  due  to  three  independent  charac- 
teristics of  the  model:  its  lift,  its  displacement,  and  the  displacement  of  its  wake.  The  contri- 
butions of  these  three  characteristics  are  usually11  estimated  by  representing  them  as  a vortex, 
a doublet  and  a source  respectively,  and  the  three  corresponding  modes  of  wall  interference 
referred  to  as  lift  interference,  solid  blockage,  and  wake  blockage.  If  the  velocity  potentials 
are  obtained  for  each  of  these  singularities  in  the  tunnel,  then  the  interference  will  be  represented 
by  interference  potentials,  which  are  the  difference  between  the  potentials  in  the  tunnel  and  in 
free  air.  Evaluation  of  these  interference  potentials  in  the  neighbourhood  of  the  model  may 
then  be  used  as  estimates  of  the  interference  corrections  required  to  be  made  to  measured 
quantities. 


3.2  Subsonic  I. incar  Theory 


The  linearized  equation  of  motion  for  two-dimensional  subsonic  compressible  flow  is 
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(3.1) 
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where  <t>  is  the  perturbation  velocity  potential  of  the  flow.  Baldwin  el  u/.15  have  derived  a single 
homogeneous  boundary  condition  for  ventilated  wind  tunnel  walls  in  the  form 


hh  I t><)> 

2 ' P iz 


% *r 


where  P is  a porosity  parameter,  and  F is  slot  parameter  given  by 
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F — In  cosec  (3.3) 
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where  u is  the  slot  width  and  d the  mean  periodic  slot  spacing  (see  Fig.  3). 

It  should  be  noted  that  the  general  boundary  condition  (3.2)  contains  the  special  cases 

(i)  solid  wall — 1 IP  > oo  or  F -*■  oo; 

(ii)  perforated  wall — \/P  ^ 0,  F = 0; 

(iii)  ideal  slotted  wall — I jP  = 0,  F ^ 0; 

(iv)  open  jet — \/P  = 0,  F 0 

and  hence  results  derived  from  this  boundary  condition  will  contain  results  applicable  to  all 
possible  tunnel  wall  conditions. 

Writing 

t = h -i  <f>*  (3.4) 

where  ^ i is  the  perturbation  potential  about  the  model  in  free  air  and  <j>*  is  the  interference 
potential,  then  Baldwin  el  a/.15  and  Wright16  have  obtained  solutions  for  <f>*  by  solving 
Equations  (3.1)  and  (3.4)  with  the  boundary  condition  (3.2)  by  a Fourier  transform  method. 
This  then  yields  the  change  in  stream  conditions  caused  by  the  walls  at  the  model  position, 
given  by 
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which  is  the  blockage  correction. 


which  is  the  upwash  correction,  and 
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which  is  the  streamline  curvature  correction.  Each  correction  will  now  be  considered  in  turn. 


3.2.1  Solid  Blockage 

In  free  air,  the  potential  due  to  a doublet  of  strength  ma  is  given  by 


and  the  strength  of  the  doublet  may  be  related  to  the  size  of  the  model  with  sufficient  accuracy  by 


where  A is  the  cnss-sectional  area  of  the  model  in  the  x-z  plane.  Baldwin  et  a!,  showed  that  the 
interference  potential  due  to  the  doublet  is  given  by 
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Differentiating  (3. 10)  and  setting  t ; Ogives 
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lor  the  additional  stream  velocity  at  the  model  position  due  to  the  walls  Allowing  h » oo 
gives  the  solid  wall  result 


3.2.2  Wake  Blockage 

For  a two-dimensional  source  of  strength  m„ 
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j 'n  (v*  | (3.14) 
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and  the  strength  of  the  source  may  he  related  to  the  model  wake  by 

m»  iUcCn  (3.15) 

where  Cn  is  the  drag  coefficient  of  the  model.  The  interference  potential  is  then  given  by 
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where  /.<  and  /<•  are  as  defined  in  (3.1 1).  Differentiating  and  setting  v ; 0 yields 
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tor  the  additional  stream  velocity  at  the  model  position  due  to  wake  blockage.  Once  again, 
allowing  P - 0 gives  the  solid  wall  result 


It  should  be  noted  that  the  expressions  (3.17)  and  (3.18)  are  different  from  those  derived 
by  Wright,1*  which  included  an  error.  This  error,  which  was  pointed  out  by  ('athcrall17  has 
been  corrected  in  the  above  equations. 
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3.2.3  I if*  Interference 

In  free  air  the  potential  ilue  to  a two-dimensional  vortex  at  the  origin  is  given  by 
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The  circulation.  I',  about  the  model  is  related  to  the  model  lift  per  unit  width  by 
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and  the  interference  potential  may  be  written  as 
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h ^sinhi/  I hq  cosh  | cosh  - q 
Differentiating  and  setting  i c 0 yields  the  upwash  velocity,  given  by 
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and  the  strcamwisc  gradient  ot  the  upwash  velocity  or  the  streamline  curvature  at  the  origin  is 
given  by 
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Again  allowing  /■' 


X'.  gives  the  solid  wall  result 
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Once  again  the  expressions  (3.23)  and  (3.24)  are  different  front  those  derived  by  Wright  due 
to  an  error  in  the  analysis  also  pointed  out  by  Cutherall. 


3.2.4  Blockage  Velocity  Gradients 

The  variation  in  longitudinal  velocity  due  to  both  solid  and  wake  blockage  given  by 
Equations  (3.12)  and  (3.17)  is  not  symmetrical  forward  and  aft  of  the  aerofoil  mid-chord.  This 
gives  rise  to  longitudinal  velocity  gradients  given  (at  the  origin)  by 


for  the  solid  blockage  term,  and 
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for  the  wake  blockage  term.  These  longitudinal  velocity  gradients  give  rise  to  longitudinal 
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bouyancy  forces  on  the  mode!  which  would  he  absent  in  free  air  Balance  measurements  of  model 
drag  would  record  this  bouyancy  force  and  would  therefore  have  to  be  corrected  in  order  to 
correspond  to  free-air  values.  However,  if  model  drag  is  obtained  by  wake  traverse  measure- 
ments downstream  of  the  model,  no  bouyancy  corrections  need  be  applied  Since  no  balance 
drag  measurements  were  undertaken  in  the  present  tests,  these  terms  will  not  be  considered 
further 


3.3  Corrections  to  Measured  Quantities 

As  shown  in  Section  3.2.3.  lift  interference  induces  an  upwash  velocity  w which  varies 
along  the  tunnel  and  necessitates  corrections  to  the  model  incidence,  loading,  and  lift  and 
pitching  moment  coefficients.  The  correction  to  incidence,  A»,  represents  the  interference  upwash 
at  the  aerofoil  mid-chord,  and  the  remaining  corrections  allow  for  the  residual  upwash  distributed 
over  the  aerofoil.  Since  the  representative  vortex  must  be  placed  at  the  aerofoil's  centre  of  pressure. 


and  the  mid-chord  is  approximately  at  a distance  < 
to  incidence  will  be 
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downstream,  then  the  correction 
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The  choice  of  Aa  is  such  that  the  residual  loading  correction.  A/.,  vanishes  at  the  aerofoil 
leading  edge,  and  the  increment  to  be  added  to  the  measured  loading  is 
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where  t 1 is  a parameter  used  to  describe  the  chordwisc  distance  aft  of  the  leading  edge  of  the 
aerofoil. 
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The  incremental  correction  to  lilt  is  obtained  by  integrating  A/  along  the  aerofoil  chord, 
giving 
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and  the  incremental  correction  to  pitching  moment  is  obtained  by  integrating 

{A/.(.vi,  | ()•  25c  v)| 

along  the  aerofoil  chord  giving 
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Ii  should  he  noted  that  the  values  of  lift  and  pitching  moment  coefficients  appearing  in 
Equations  ( f 30),  (3  34)  and  ( t 33)  arc  uncorrected  and  based  on  the  uncorrected  kinetic  pressure, 
kpl* 

The  doublet  and  source  representing  the  aerofoil  displacement  and  that  of  Us  wake  give 
rise  to  changes  in  longitudinal  velocity  due  to  solid  and  wake  blockage  respectively  these  two 
blockage  components  can  be  added  together  to  give  a total  blockage  factor 


or 


where 


It  is  usual  to  express  the  values  of «,  and  <„  for  the  general  ventilated  wall  case  by  factoring 
the  solid-wall  blockage  factors  by  the  quantities  t),  and  W„  respectively  I he  total  blockage 
factor  can  then  be  written 
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where  *»,■  and  «Mr  are  the  solid  wall  values  of  t,  and  •„  and  W,  and  tl„  are  the  ratios  of  solid  wall 
to  general  ventilated  wall  blockage  factors.  Thus  11,  ll„  I for  the  solid  wall  case  I rom 
Equation  (3  13), 
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and  from  Equation  (3  IK). 
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where,  if  drag  measurements  are  made  by  the  wake  traverse  method,  Co  is  the  uncorrected 
drag  coefficient  based  on  the  uncorrected  kinetic  pressure  jr»C2.  Garner  «•/  <r/."  point  out  that 
improved  agreement  with  experimental  results  is  obtained  if  «,r  and  «„,•  are  modified  on  the  basis 
of  more  elaborate  theories  Thus  «„■  contains  a factor  to  account  lor  the  effect  of  incidence, 
and  the  cross-sectional  area  is  multiplied  by  an  empirical  factor. 


(3  41) 


The  expression  for  ««,<■  contains  an  extra  compressibility  factor  which  then  gives  a com- 
pressibility factor  between  (f  - and  ft  :l,  thus 
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and  this  expression  has  been  found  to  lead  to  better  agreement  with  experiment  at  high  sub- 
sonic Mach  numbers. 


Because  blockage  changes  the  stream  longitudinal  velocity,  the  static  pressure,  density 
and  static  temperature  of  the  stream  will  also  change.  The  incremental  corrections  to  be  added 
to  the  undisturbed  stream  values  to  allow  for  these  changes,  and  the  associated  corrections  to 
flow  quantities  derived  from  them  are  as  follows: 
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It  is  now  possible  to  correct  aerixlvnumic  ci>ellicieiils  measured  in  tlie  tunnel,  tine1  based 
on  uncoiieited  kinetic  piessuie,  to  lice  an  values  line  ptixrss  can  be  considered  ill  two  steps 
(il  adding  ti>  the  nicuMiied  Ci  and  ( values  the  lit  t interference  correct  nuts  ACi  and 
-\<  . and 

tut  multiply  ill);  the  results  bv  (»'.  the  ralio  ol  uncorrected  to  corrected  kinetic  pressure 
1 he  lice  an  values  ol  aerixly  naniic  civellicients,  based  on  the  corrected  kinetic  pressure  are 
then  given  bv 
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When  comparing  results  (tout  tests  on  two  01  more  models  01  liom  tests  m dittoing  wall  con 
tigurations.  corrected  values  ol  both  Mach  number  and  angle  ot  incidence  must  agree  lot  the 
comparison  to  be  valid  In  practice,  tins  seveielv  limits  the  inimhei  ol  cases  available  lot  coni 
panson  One  method  ol  overcoming  this  limitation  is  to  use  zero  lilt  lilt  curve  slo|x-  as  the  basis 
for  companson.  in  which  case  the  onlv  lecpniement  is  that  comparisons  lx-  made  at  the  same 
corrected  Mach  niimbei  It  is  therefore  useful  to  combine  the  separate  corrections  to  lift 
cix'thcieiu  and  angle  ol  incidence  into  a single  correction  to  lilt  curve  slo|x-  Assuming  that,  in 
a small  region  around  zero  lilt,  the  lift  curve  is  approximately  lineal,  i e 
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then  the  lift  curve  slo|x'  mav  lx-  written  as 
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and  the  lilt  curve  slo|x-  corrected  to  fteeair  conditions  as 
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Substituting  into  ( t s ')  we  obtain 
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From  (3.30)  wc  have 
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Assuming  lhai  the  pitching  moment  coefficient  varies  linearly  with  incidence  in  the  region  of  the 
zero  lift  point,  i.e. 
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Also,  from  (1 .14) 
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Substituting  (3  56)  and  (3  57)  into  (3.54)  we  obtain 
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I his  expression  takes  into  account  both  blockage  and  lilt  interference  and  thus  enables  lift  curve 
slopes  to  be  directly  compared,  the  only  further  requirement  being  that  the  individual  lift  curve 
slopes  be  compared  at  their  corrected  Mach  numbers. 


3.4  (General  Remarks 

I he  theory  of  wind  tunnel  wall  interlerence  correction  is  still  far  from  complete.  The  basic 
deficiency  in  the  application  of  subsonic  theory  to  ventilated  wall  test  sections  is  an  inadequate 
knowledge  ol  the  wall  boundary  conditions  and  a lack  of  understanding  of  the  viscous  loss 
characteristics  across  the  walls.  I his  makes  it  difficult  to  apply  classical  wall  correction  theory 
and  expect  the  same  accuracy  that  has  been  realized  Tor  solid  walls  or  for  open  |cts  If  the  wall 
boundary  conditions  lor  ventilated  test  sections  could  be  specified  as  satisfactorily  as  those  for 
solid  walls  or  lor  open  jets,  ventilated  wall  corrections  would  be  subject  to  the  same  limitations 
that  have  always  been  present  in  the  classical  correction  methods  These  limitations  arise  from  the 
lirst  order  linearized  approach  to  the  problem  and  to  the  modelling  of  the  situation  by  three 
independent  corrections  This  gives  rise  to  the  implication,  for  example,  that  since  there  is  stream- 
line curvature,  a warped  model  should  be  tested  to  simulate  the  uniform  unbounded  flow  situation, 
and  that  the  model  warpage  should  be  changed  with  each  angle  of  attack.  However,  the  wall  loss 
characteristic  problem  is  more  fundamental  to  the  ventilated  wall  case  since  it  precludes  the 
confident  application  of  classical  methods. 

C onsider  the  slot  parameter,  F,  defined  by  Equation  (3.3).  It  was  pointed  out  by  Davis  and 
Moore,1®  and  later  by  Bcrndt  and  Sftrensen,'2  that  if  the  slots  have  finite  thickness,  / (see  Fig. 3). 
(hen  F is  given  by 
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More  recently  Barnwell-"  has  generalized  the  flat  slat  formulation  so  that  variations  involving 
slots  with  sidewall'  and  separation  in  the  plenum  can  he  treated  Me  derived  an  expression  for  A 
for  the  case  in  which  the  flow  separated  on  the  plenum  side  of  the  slots  as  follows: 


and  lor  the  case  of  slots  with  sidewalls  (see  Fig  31 
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It  can  he  shown  that  the  differences  in  the  values  of  A derived  from  l quations  (3.3).  (3  60)  or 
(3.61)  arc  relatively  small  if  the  slot  width  to  slot  spacing  ratio  a ,1  is  small.  However,  these 
differences,  due  to  the  uncertainty  of  the  exact  nature  of  the  wall  boundary  condition,  can 
become  quite  large  in  the  case  of  a tunnel  wall  with  only  a few  relatively  wide  slots 

In  Reference  21  ( hen  and  Mears  develop  a model  in  which  the  slots  are  represented  by 
bodies  which  arc  generated  by  doublet  rods  of  constant  strength  as  shown  in  figure  4 The 
boundary  of  each  knly  is  the  dividing  streamline  between  the  outside  (low  and  the  (low 
associated  with  the  doublet  rod  lhe  slot  width  and  slot  spacing  are  again  designated  as  u and  J. 
and  the  distance  between  the  doublet  rinls  is  designated  as  />  Barnwell-"  has  pointed  out  a 
mistake  in  the  analysis  of  Chen  and  Mears  in  that  they  equated  the  distance  between  the  rods.  /). 
with  the  slot  width  a.  These  distances  are  obviously  not  equal  since  there  must  be  a small  space 
between  the  end  of  the  rod  and  the  dividing  streamline  as  indicated  in  figure  4u  so  that  the 
fluid  which  is  emitted  in  front  of  the  ri>d  by  the  source  half  of  each  doublet  can  flow  around  the 
end  and  into  the  sink  half  of  the  doublet  at  the  rear.  This  error  led  to  an  order  of  magnitude 
difference  in  the  values  of  A for  small  values  of  a J and  to  some  doubt  as  to  the  validity  of  the 
approach  when  compared  with  the  classical  approach  leading  to  (3  3). 

Chen  and  Mears  analysis  leads  to  an  expression  for  A'  of 
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where  in  this  case  / is  the  centreline  thickness  of  the  slat.  The  rod  gap  />  is  unknown  and  must 
be  determines!  iteratively  by  fixing  J and  / and  assuming  a value  for  P.  The  value  of  a is  then 
obtained  from  the  equation  for  the  stream  function  describing  the  slat  shape: 
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where  /i  is  related  to  the  doublet  strength  and  is  given  by 
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The  slat  shapes  developed  by  Chen  and  Mears-  approach  are  plotted  in  figure  4/>  for 
several  values  of  centre-line  thickness  It  is  apparent  that  the  effective  thickness  of  the  slat  cannot 
be  reduced  below  about  J.  For  this  reason,  the  usefulness  of  the  doublet  rod  model  in  its  present 
form  is  limited.  More  control  could  be  exerted  over  the  effective  thickness  of  the  slat  if  the 
doublet  strength  could  be  varied  along  the  length  of  the  rod  However,  the  model  does  indicate*" 
that  the  value  of  A’  is  more  sensitive  to  the  shape  of  (he  slat  than  to  whether  or  not  the  flow 
separates  on  the  plenum  side  of  the  slat. 

A more  serious  problem  in  the  specification  of  the  wall  boundary  condition  arises  with  the 
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term  to  take  account  of  viscous  effect'*  in  the  slots  The  constant  "porosity",  /’.  in  Fquution  (3.2) 
is  assumed  to  he  determined  by  the  wall  geometry  and  the  I'ree-streatn  Mach  number,  although 
unlike  the  slot  parameter,  / . there  is  no  method  for  computing  it  a priori.  This  term  is  analogous 
to  the  term  expressing  the  linearized  pressure  drop  across  a porous  wall  Ihe  basis  for  using 
such  a characteristic  for  slotted  walls  is  however,  even  more  tenuous  than  the  basis  for  using  it 
lor  porous  walls,  slotted  walls  departing  even  further  from  the  physical  basis  (i.e  Darcy's  law 
for  flow  through  dense  porous  media).  It  has  been  suggested  by  several  workers1-  that  the 
viscous  loss  characteristic  of  slotted  walls  may  be  better  represented  by  a quadratic  term  rather 
than  a linear  one.  especially  in  cases  where  the  cross-flow  velocities  are  large 

I he  whole  area  of  slotted  wall  boundary  conditions  is  further  complicated  by  the  presence 
of  the  wall  boundary  layers,  no  account  of  which  is  taken  by  representations  such  as  (3.2)  The 
total  cross-sectional  area  of  the  wall  boundary-layers  may  well  approach  or  even  exceed  the 
cross-sectional  area  of  the  model  It  is  therefore  possible  that  to  a large  extent  the  fluid  entering 
the  slots  could  lx-  low  momentum  boundary-layer  air,  rather  than  free-stream  air  The  slot  flow 
is  then  bound  to  he  less  efficient  in  maintaining  a pressure  difference  between  the  test  section 
and  the  plenum  chamber  1 hese  considerations  also  apply  to  any  region  where  the  slot  flow 
passes  into  the  test  section  from  the  plenum  chamber  with  its  more  or  less  stagnant  air  l or  the 
case  where  there  is  vanishing  longitudinal  flux  in  the  slot  region,  there  will  be  a negligible  pressure 
difference  and  the  wall  will  appear  to  the  test  section  flow  as  a free-jet  boundary  at  the  plenum 
pressure 

It  is  clear  then,  that  the  flow  properties  of  slotted  walls  may  well  be  rather  complicated,  very 
likely  more  complicated  than  consistent  with  the  simplifying  assumptions  made  in  the  homo- 
geneous boundary  conditions  of  classical  interference  theory  presented  in  the  previous  section. 
Some  attempts'-  have  been  made  to  examine  the  situation  with  a view  to  obtaining  a better 
representation  of  the  boundary  conditions,  but  these  await  further  experimental  verification 
before  the  use  of  then  more  complicated  forms  can  be  justified.  Meanwhile  the  use  of  the  classical 
forms  of  boundary  condition  will  continue,  but  the  accuracy  of  the  resulting  corrections  will  be 
open  to  doubt  and  will  certainly  be  less  than  that  achieved  for  the  solid  wall  or  open  jet  cases. 

4.  APPLICATION  Ol  UM  AR  TIIFOR3 
4.1  Solid  Wall  Results 

I he  results  of  tests  on  both  models  of  NACA-0012  and  BCik- 1 in  the  solid  wall  test  section 
have  been  analysed  in  accordance  with  linear  interference  theory  presented  in  the  previous  section. 

I ift  curve  slope  at  zero  lift  was  used  as  a basis  for  comparison,  the  data  being  corrected  by  the 
application  of  I quation  (3  58)  I or  the  solid-wall  case,  the  lift  interference  parameters  iS„  and 
<S i are  given  explicitly  by  ( 3 3 1 ).  with  Sn  going  to  zero  The  blockage  factors  «„■  and  * «T  are  given 
by  (3  41)  and  (3  43)  No  wake  drag  measurements  were  undertaken  for  the  solid  wall  tests,  and 
the  valir  of  drag  coefficient  in  (3  42)  was  therefore  obtained  from  the  results  of  the  slotted  wall 
tests  at  the  corresponding  uncorrected  Mach  numbers  and  attitudes.  Since  wake  blockage 
accounts  for  less  than  10  ,,  of  the  total  blockage  factor,  this  should  not  give  rise  to  any  significant 
errors. 

I he  corrected  values  of  lilt  curve  slope  are  plotted  in  figure  5 for  BGK-1  and  in  Figure  b 
for  the  NACA-0012.  Results  from  tests  on  the  two  models  of  each  section  shape  now  show 
satisfactory  agreement,  although  scatter  is  still  evident  in  the  results  for  the  BC>k-l  section, 
especially  at  supercritical  Mach  numbers.  The  corrected  data  generally  fall  within  ; 2",,  of  the 
common  curve,  compared  with  differences  of  about  I0"„  for  the  uncorrected  data  (see  Figs. 
7 and  S)  An  analysis  of  the  expected  errors  in  the  value  of  lift  curve  slope  (on  the  basis  of  typical 
experimental  errors  in  angle  of  incidence  of  0 02  and  in  lift  coefficients  of  0 002,  and  lift 
curve  slopes  being  derived  from  the  slope  of  the  lift  curve  over  an  angle  of  incidence  of  4 ) 
indicates  that  they  would  lx-  of  the  order  of  - 2"„  the  same  order  as  the  scatter  between 
corrected  results. 

f urther  evidence  of  the  success  of  linear  theory  in  correcting  the  data  obtained  in  the  solid 
wall  tests  is  provided  hv  the  zero  lift  NACA-0012  results.  Figure  presents  the  variation  of 
pressure  ratio  with  Mach  number  for  one  particular  chordwise  position  ( \ c 0-3).  Corrected 
data  from  both  models  of  the  NACA-0012  are  plotted  and  the  data  are  seen  to  define  a single 
curve  with  very  little  scatter. 
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FIG.  6 VARIATIONS  OF  LIFT  CURVE  SLOPE  WITH  MACH  NUMBER  - NACA  0012 

CORRECTED  SOLID  WALL  DATA 


FIG.  7 VARIATION  OF  LIFT  CURVE  SLOPE  WITH  MACH  NUMBE 

UNCORRECTED  SOLID  WALL  DATA 


UNCORRECTED  SOLID  WALL  DATA 


FIG.  9 VARIATION  OF  PRESSURE  RATIO  WITH  MACH  NUMBER  - NACA  001? 

CORRECTED  SOLID  WALL  RESULTS 
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As  suggested  in  Relerencc  the  existence  of  a design  pressure  distribution  tot  the  Itiik  1 
section  also  pros  isles  a conxenient  reference  point  tor  the  studs  of  ssmd  tunnel  wall  interference 
I lie  experimental  design  distributions  base  been  found  to  he  quite  sensitise  to  small  changes 
in  W and  t.  and  such  a reference  point  would  appear  to  be  well  defined  It  the  results  could  be 
nude  tree  ol  interference  elfects,  the  conditions  at  which  the  design  pressure  distribution  was 
obtained  for  each  model  should  he  the  same  (since  viscous  effects  should  be  in  the  same  order 
for  h’th)  although  this  condition  would  not  he  expected  to  agree  with  that  found  from  insiscid 
tlieois  Applsing  linear  thcorv  to  the  lUik  I design  pressure  distribution  results  from  the  solid 
wall  tests  we  obtain  the  collected  design  point  conditions  as  tabulated  below 
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Again  the  agreement  between  the  two  models  is  excellent  I he  existence  of  a single  design 
condition  relevant  to  both  models  is  further  supported  bv  the  close  agreement  between  the 
experimental  design  pressure  distributions,  as  shown  in  I igure  lit 

Overall,  the  solid  wall  data  corrected  In  linear  tlieorv  shows  \et\  salislactorv  agreement 
It  is  theretorc  to  be  expected  that  the  cot  reeled  results  will  pros  ule  Ian  Iv  close  to  interference  free 
data 

4.2  Molted  Wall  Results 

In  the  slotted  wall  case,  the  tunnel  walls  are  characterized  In  two  non-dimensional 
parameters  / which  is  a unique  function  of  the  geometrv  of  the  slots,  and  /’  (the  “porositv 
lactor  1 usuallv  considered  as  the  ratio  ri  /’.  which  appears  in  the  full  homogeneous  bouintarv 
condition  (52)  to  allow  lor  v iscous  effects  in  the  slots  1 he  parameter  I is  faith  easih  determined, 
the  onh  problem  being  whether  to  take  the  timte  thickness  of  the  slats  mto  account  /'  01  S /’ 
.m  the  other  hand  must  be  determined  empiricallv  and  us  definition  is  the  basic  problem 
encountered  in  the  application  ot  subsonic  linear  tlieorv  to  slotted  wall  wind  tunnels 
l otisidering  tirstlv  the  slot  geometrv  parameter  / . from  I quatton  |i  t| 


l or  the  slotted  test  section  walls  used  in  the  present  tests 

ii  12  7 mm 
b 55.1  mm 
h SOS  m ni 

and 


<1  b \ ’<>  I mm  ( \ ’1 

gi'  mg 

/ 0 OS  10 

It  the  finite  thickness  ol  the  slats  (/  7 54  111111I  is  taken  into  account,  then  from  Imitation 

1 1 5d)  we  obtain  I 0 1 Id 

I timing  now  to  the  \ iscous  How  parameter  /\  there  arc  three  classes  of  methods,  all 
basualK  experimental,  which  haxe  boon  used  t v>r  its  determination  tho\  are 
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FIG.  10  COMPARISON  OF  BGK  1 DESIGN  PRESSURE  DISTRIBUTIONS; 

SOLID  WALL  DATA 


direct  measurement; 
indirect  measurement . 
comparison  of  measurements. 

Direct  measurement  consists  of  the  calibration  of  actual  or  sample  wind  tunnel  walls  by  the 
measurement  of  the  pressure  drop  across  the  wall  for  various  cross  llow  velocities  and  free 
stream  Mach  numbers.  I his  method  is  more  suitable  lor  use  with  perforated  walls,  and  has 
been  extensively  used  for  this  purpose  at  the  Arnold  fngineering  Development  Centre  (AFDC).'-2 
A method  of  indirect  measurement,  originating  at  the  National  Lueht-en  Ruimtevaart- 
laboratorium  (NIR),-'J  consists  of  an  aerofoil  drag  balance  method  in  which  all  non- 
interference drag  contributions  are  either  measured  or  calculated,  and  are  then  compared  to 
the  drag  obtained  from  wake  rake  measurements.  The  "porosity  factor"  ft  /’  may  then  be 
obtained  for  each  Mach  number  and  angle  of  attack  by  subtraction.  This  method  sutlers  from 
the  disadvantages  of  requiring  unrealistically  accurate  measurements  of  profile  drag  to  be 
obtained  from  aerofoil  surface  pressure  measurements  and  the  difficulty  of  accurately  estimating 
skin-friction  drag 

The  methods  allowing  the  determination  of  ft  /’  from  comparative  measurements  usually 
consist  of  either  testing  a single  aerofoil  model  in  a test  section  in  which  the  wall  conditions 
may  be  varied  front  closed,  through  several  wall  permeabilities,  to  a condition  comparable  to 
that  of  an  open  jet,  or  of  testing  two  or  more  geometrically  similar  models  of  the  same  aerofoil 
section  in  a single  ventilated-wall  test  section.  Such  methods,  because  of  their  straightforward 
application,  are  by  far  the  most  popular. 

I he  approach  adopted  here  will  also  be  one  of  comparative  measurements,  hut  because 
of  the  unusual  set  of  test  results  available,  the  approach  will  be  somewhat  different  from  that 
usually  followed  In  the  previous  section  it  was  shown  that  the  solid  wall  results,  corrected  accord- 
ing to  linear  interference  theory,  provide  a set  of  measurements  close  to  being  free  of  the  effects 
of  interference.  The  approach  will  be  to  accept  these  measurements  as  being  truly  "interference 
free"  and  to  use  them  to,  in  effect,  calibrate  the  slotted  walls  in  finding  the  applicable  value  of 
fii'P.  The  basis  for  comparison  will  be  the  variation  of  corrected  solid  wall  lilt  curve  slope  (at 
zero  lift)  with  Mach  number,  and  then  to  attempt  to  find  a value  of  /’  which  minimizes  the 
deviations  of  the  corrected  slotted  wall  lilt  curve  slopes  from  these  "interference  free"  measure- 
ments for  this  purpose  "mean  curves"  shown  on  figures  5 and  6 have  been  drawn  as  lines 
of  best  lit  to  the  corrected  solid  wall  data.  To  quantify  the  agreement  between  the  slotted  wall 
data  and  these  mean  curves,  we  define  a root  mean  square  deviation  (A™,*)  given  by 


where  K,  t are  Ihe  corrected  solid  wall  lift  curve  slopes  taken  from  the  mean  curve,  and  A,  are 
the  slotted  wall  lilt  curve  slopes,  the  summation  being  taken  over  the  n values  of  Mach  number 
at  which  measurements  were  made  for  each  model  of  each  aerofoil  section.  Care  was  taken  to 
ensure  that  the  value  of  A',/  was  interpolated  from  the  mean  curves  to  correspond  to  the  corrected 
Mach  number  of  the  slotted  wall  values  AY.  I ach  model  of  each  aerofoil  section  was  treated 
separately  to  ensure  that  any  trends  due  to  either  varying  chord  or  aerofoil  section  would  be 
apparent. 

As  for  the  solid  wall  results,  the  lift  curve  slopes  were  corrected  by  the  application  of 
liquation  (3.58).  In  this  case,  however,  the  interference  parameters  <S()  and  iq  given  by  (3.31) 
now  become  functions  of  /■'  and  ft  /’  and  are  found  by  substituting  from  (3.23)  and  (3.24)  into 
(3.31).  Similarly  the  blockage  factor, ««,  w hich  determines  (l  (ihe  ratio  of  corrected  to  uncorrected 
kinetic  pressure),  is  related  to  /•'  and  ft  /’  via  liquations  (3  12).  (3. 1 7)  and  (3.37). 

The  variation  of  the  root  mean  square  deviation  Arms  with  ft  /’  is  shown  in  f igure  II  for 
the  four  aerofoil  models.  I he  values  of  ft  /’  giving  the  minimum  deviation  for  each  model  are 
also  noted  on  the  figure.  Although  there  is  some  variation  in  this  value  between  the  models, 
notably  for  the  101  6 mm  chord  NACA-0012  model,  the  results  indicate  a single  value  of  ft  /’. 
appropriate  to  the  slotted  walls  used,  with  a value  of  about  0-4  It  is  interesting  to  note  the 
steeper  slopes  about  the  minima  obtained  for  both  the  larger  chord  models.  This  is  indicative 
of  the  greater  dependence  of  the  larger  chord  model  results  on  fi  /’  due  to  the  presence  of  increased 
i n t e r fe  re  nee  e ffee t s . 


These  calculations  were  repeated  with  the  value  of  / (0  II1))  taking  into  account  the  finite 
thickness  of  the  sluts.  I he  indicated  value  of ' /’  differed  by  only  a few  per  cent  from  the  original 

value  (0  4)  and  the  value  of  /•  was  therefore  retained  at  0 081  for  all  future  work. 

Corrected  lift  curve  slopes  were  then  calculated  using  linear  theory  with  the  value  of  ff  /' 
held  constant  at  0-40  These  corrected  values  arc  plotted  in  Figure  12  for  the  BciK-l  section 
and  in  Figure  If  for  the  NACA-0012  together  with  the  mean  curves  of  the  corrected  solid  wall 
results  for  comparison.  Apart  from  the  low  Mach  number  range  for  the  BtiK-1  section,  the 
slotted  wall  results  now  tall  within  • 2 to  3"„  of  the  corrected  solid  wall  curves.  Considering 
that  the  increase  in  lilt  curve  slope  over  the  uncorrected  results  (which  are  shown  in  Figures 
14  and  15  for  comparison)  amounts  to  between  25“,,  and  50",,  for  the  203 '2  mm  chord  models 
and  between  15",,  and  25",,  lor  the  101  b mm  chord  models,  this  agreement  is  remarkably  good. 

Once  again,  the  zero  lilt  NACA-0012  results  provide  a useful  check  on  the  blockage  correc- 
tions in  the  form  of  pressure  ratio  (/>  //)  variation  with  Mach  number  at  constant  chordwise 
position.  If  we  again  take  the  corrected  solid  wall  data  as  being  interference  free,  and  provided 
that  a given  pressure  at  a fixed  point  on  the  model  represents  a given  flow  on  the  model,  i.e. 
provided  the  blockage  effect  produces  no  distortion  at  the  model  and  can  he  applied  as  a correc- 
tion to  longitudinal  velocity,  then  the  differences.  V)/,  on  horizontal  intercepts  between  the 
corrected  solid  wall  curve  and  the  uncorrected  siotted  wall  curve  give  the  blockage  effects  on 
Mach  number  for  the  slotted  walls.  A typical  plot,  for  v c 0 03  is  shown  in  Figure  lb.  From 
this  plot,  at  a Mach  number  of  0-75  we  sec  that  SM  is  approximately  0-009,  compared  with 
0 002  predicted  by  linear  theory,  linear  theory  clearly  underestimates  the  blockage  effect 
by  a significant  amount,  and  this  is  true  throughout  the  Mach  number  range. 

Hy  plotting  these  A l/s  against  Mach  number  we  may  obtain  an  estimate  of  the  actual 
blockage  effects  on  Mach  number  as  a function  of  Mach  number  for  the  slotted  walls.  If  similar 
curves  are  plotted  for  various  chordwise  positions  a check  can  be  made  for  the  existence  of 
significant  flow  distortions,  since  in  that  case  it  would  be  expected  that  the  shape  of  these  curves 
would  vary  in  an  ordered  manner  with  chordwise  position.  Such  a series  of  results  are  plotted 
in  Figure  17,  and  although  significant  scatter  is  present,  the  scatter  is  largely  random,  the  result 
of  experimental  errors  rather  than  any  trend  due  to  significant  flow  distortion,  supporting  the 
assumption  that  blockage  can  be  allowed  for  by  a simple  correction  to  tree-stream  velocity. 
Also  plotted  on  the  figure  arc  the  curves  predicted  by  linear  theory  for  the  case  of  a slotted  wall 
with  fi  r 0 40  and  an  open  jet.  It  is  evident  that  even  the  open  jet  case  underestimates  the 
blockage  hy  a factor  of  four  or  five.  This  would  suggest  that  the  present  slotted  walls,  at  least 
as  far  as  blockage  is  concerned,  are  more  "open"  than  an  open  jet.  Similar  results  have  been 
found  by  Penrcey  <-/  at. - 1 and  some  possible  reasons  for  this  discrepancy  will  be  discussed  in  a 
later  section. 

Classical  linear  theory  suggests  that  the  variation  with  Mach  number  of  the  differences 
shown  in  Figure  17  should  be  of  the  form 
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where  <■  is  an  empirical  blockage  factor.  This  empirical  blockage  factor  therefore  replaces  < n in 
Equation  (3.45)  and  embraces  the  solid  blockage  and  wake  blockage,  where  it  exists,  both  now- 
factored  by  l//J:l.  Again  following  classical  theory,  we  will  assume  that  <•  is  related  to  the  dimen- 
sions of  the  test  section  and  model  by 
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where  -f,-  is  the  effective  model  cross-sectional  area  given  by  (see  Equation  (3,41)) 
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and  F should  be  a constant  depending  only  on  the  properties  of  the  particular  test-section  walls. 
Equation  (4.2)  is  plotted  in  Figure  17,  and  the  value  of  /.'  giving  the  closest  agreement  with  the 
data  is  about  I -25.  Unfortunately,  the  data  from  the  101  -6  mm  chord  model  of  the  NACA- 
0012.  which  could  have  provided  an  independent  check  on  the  value  of  / from  a similar  process, 
contained  too  much  experimental  scatter  to  allow  any  meaningful  deductions.  (It  should  be 
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FIG.  12  VACATION  OF  LIFT  CURVE  SLOPE  WITH  MACH  NUMBE 

CORRECTED  SLOTTED  WALL  DATA 


CORRECTED  SLOTTED  WALL  DATA 


FIG.  15  VARIATION  OF  LIFT  CURVE  SLOPE  WITH  MACH  NUMBER  - NACA  0012; 

UNCORRECTED  SLOTTED  WALL  DATA 


FIG.  16  VARIATION  OF  PRESSURE  RATIO  WITH  MACH  NUMBER  - NACA  0012; 

COMPARISON  OF  SLOTTED  WALL  DATA  WITH  CORRECTED  SOLID  WALL  DATA 
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FIG.  17  VARIATION  OF  MACH  NUMBER  INCREMENT  WITH  MACH  NUMBER  - NACA  0012' 

SLOTTED  WALL  DATA 


noted  that  the  A Ms  predicted  from  Equation  t4  2)  would  he  a factor  of  4 smaller  than  those  of 
the  203-2  nun  chord  model  ) It  is  worth  noting  that  values  of  A corresponding  to  solid  walls 
and  to  slotted  walls  with  a value  ot  /’ot  0 4 (calculated  by  assuming  the  value  of  wake  blockage 
relevant  to  a Mach  number  ol  0 '51  are  0 50  and  O f respectively  and  for  the  case  of  an 
open  jet  (where  the  wake  blockage  is  zero)  the  value  of  A'  is  0 26  It  is  seen  that  the  empirical 
value  of  A.  I 25.  is  almost  5 times  the  open  iet  value 

This  empirical  blockage  correction,  in  the  form  of  Equation  (4  2).  mat  be  combined  with 
the  normal  linear  theory  for  lift  interference  to  form  a "semi-empirical"  correction  scheme  for 
the  slotted  wall  results  1 he  results  of  the  application  of  this  scheme  to  correct  the  zero-lilt  lift 
curve  slopes  are  shown  in  figure  IS  where  the  root  mean  square  deviations  from  the  means 
ot  the  corrected  solid  wall  results  are  plotted  against  d /’  Although  there  is  somewhat  more 
scatter  in  the  values  ot  ' I'  giving  minimum  deviation,  plots  of  the  corrected  lift  curve  slopes 
shown  in  I igures  l’>  and  20.  using  an  average  value  of  a /’  of  0 51.  show  similar  scatter  to  those 
corrected  by  complete  linear  theory 

As  tor  the  solid  wall  results,  the  experimentally  determined  design  pressure  distributions 
lound  lor  the  Bills  I section  may  again  be  used  as  a check  on  the  validity  of  the  corrections. 
The  design  point  data,  obtained  from  Reference  7 have  been  corrected  by  the  application  of  the 
•‘semi-empirical”  scheme  developed  above,  with  A 1 - 25  and  .3  P 0 51  (The  values  of 
lilt  interference  parameters  6,,  and  i>i  obtained  with  such  a scheme  are  016  and  0 05 
respectively.)  I he  corrected  slotted  wall  results,  together  with  the  corrected  solid  wall  results 
for  comparison,  are  tabulated  below 
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0-113 

203  2 

0 748 

o 55 

0 560 

0- 124 

0-755  0 -56  ; 0-  *36 

0 118 

.Again  the  agreement  between  the  two  models  is  quite  good,  as  is  the  agreement  between 
the  corrected  results  from  both  solid  and  slotted  wall  tests  The  agreement  between  the  pressure 
distributions  for  the  experimental  design  condition  of  each  model,  plotted  in  Figure  21,  is  not  as 
good  as  that  between  the  two  solid  wall  distributions  Evidence  of  a loss  of  lift  from  a decrease 
in  lower  surface  pressure  is  apparent  for  the  203-2  mm  chord  model  up  to  mid  chord,  with  an 
increase  in  lift  from  higher  upper  surface  suction  apparent  thereafter.  This  behav  iour  is  consistent 
with  the  presence  of  significant  streamline  curvature  for  this  model.  It  is  surprising  that  such  an 
effect  should  be  visible  in  the  slotted  wall  results  and  not  in  those  from  the  solid  walls,  since  the 
relative  magnitudes  of  the  lilt  interference  parameter  Si,  which  characterizes  streamline  curvature, 
would  indicate  that  the  streamline  curvature  in  solid  walls  should  be  more  than  double  that  in 
the  slotted  walls  This  could  be  interpreted  to  mean  that  the  streamline  curvature  in  the  slotted 
walls  is  seriously  underestimated  by  linear  theory,  but  the  amount  of  data  presently  available 
does  not  allow  any  firm  conclusion  to  be  made.  To  ensure  that  this  effect  was  not  limited  to  the 
results  from  the  BGK-I  section,  a similar  comparison  was  sought  from  the  VAC.A-0012  section 
results  Such  a comparison  is  complicated  by  the  necessity  of  finding  two  sets  of  results  with  the 
same  corrected  values  of  Mach  number  and  incidence.  Two  such  cases  (for  \ft  0 -’67),  are 
plotted  in  Figure  22.  and  although  not  to  the  same  extent  as  in  the  Blik-I  case,  evidence  of 
streamline  curvature  is  indeed  present,  masked  somewhat  bv  the  slightly  higher  corrected  incidence 
and  hence  slightly  higher  lift  for  the  slotted  wall  result  1 he  agreement  in  corrected  Mach  number 
for  these  two  cases  also  provides  an  opportunity  for  the  comparison  of  the  corrected  variation  of 
lift  coefficient  with  incidence  I igure  23  shows  quite  good  agreement  between  the  results  from 
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NACA  0012  101.6  mm  Chord  0.39 
N ACA  0012  203.2  mm  Chord  0.52 
BGK  1 101.6  mm  Chord  0.53 

BGK-1  203.2  mm  Chord  0.60 


FIG.  18  VARIATION  OF  ROOT  MEAN  SQUARE  DEVIATION  WITH  0/P' 
LINEAR  THEORY  WITH  EMPIRICAL  BLOCKAGE 
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FIG.  19  VARIATION  OF  LIFT  CURVE  SLOPE  WITH  MACH  NUMBER  - BGK  1 ; 

SLOTTED  WALL  DATA  WITH  EMPIRICAL  BLOCKAGE  CORRECTION 
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21  COMPARISON  OF  BGK-1  DESIGN  PRESSURE  DISTRIBUTION- 

SLOTTED  WALL  DATA 
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FIG.  22  COMPARISON  OF  CORRECTED  PRESSURE  DISTRIBUTIONS 

NACA  0012;  M = 0.767 


solid  and  slotted  wall  tests  1 he  uncorrected  curves  are  also  included  for  comparison.  Several 
similar  opportunities  for  comparison  ol  corrected  lift  curves  are  available  front  the  BCiK-l 
data.  Unfortunately  no  cases  evist  where  the  corrected  Mach  numbers  agree  between  tests  on 
models  of  different  chords  C orrected  lilt  curves  for  those  cases  where  corrected  Mach  numbers 
are  closely  in  agreement  are  plotted  in  f igure  24  Again  the  agreement  is  quite  good  fortunately 
each  case  provides  one  pair  of  results  for  which  the  corrected  angles  of  incidence  are  very  nearly 
identical  The  pressure  distributions  corresponding  to  these  pairs  are  compared  in  Figures  25, 
2(i  and  27.  I he  agreement  between  the  pressure  distributions  with  a corrected  Mach  number  of 
0-705  (Fig.  25)  is  quite  remarkable,  the  shape  and  extent  of  the  short  supersonic  region  being 
reproduced  very  accurately.  The  comparison  of  Figure  26  is  interesting,  since  the  conditions  of 
the  comparison  a Mach  number  just  below  the  experimentally  determined  design  value  and 
an  angle  of  incidence  just  above  are  in  a region  in  which  the  shape  of  the  pressure  distributions 
are  found  to  be  extremely  sensitive  to  small  changes  in  Mach  number  and  angle  of  incidence.7-® 
The  agreement  obtained  in  this  case  significantly  increases  confidence  in  the  belief  that  the 
corrections  are  producing  truly  interference  free  results.  1 he  final  pair  ( Fig.  27),  for  a significantly 
supercritical  Mach  number  (0-765).  again  show  excellent  agreement  (considering  the  small 
difference  in  angle  of  incidence).  The  fact  that  classical  interference  theory,  which  is  not  valid 
under  supercritical  conditions,  still  produces  such  good  agreement  in  this  case  is,  perhaps,  one 
of  the  more  surprising  results  of  the  investigation 

As  a final  check  on  the  v alidity  of  the  interference  corrections,  the  chordvvise  position  of  the 
upper  surface  shock  on  the  NACA-0012  at  zero  incidence  may  lx-  used.  Figure  25  shows  the 
shock  positions  (arbitrarily  defined  as  the  mid-point  of  the  surface  pressure  rise)  for  both 
models  of  the  NACA-0012  in  both  solid  and  slotted  walls,  plotted  against  uncorrected  Mach 
number  Although  there  is  significant  scatter  present,  there  are  obvious  trends  towards  increasing 
deviations  with  increases  in  both  model  size  and  Mach  number.  In  Figure  24  the  shock  locations 
are  rcplotted  against  Mach  number  corrected  by  the  application  of  classical  theory  to  the  solid 
wall  results,  and  Equations  (4.2)  and  (4.7)  with  / l - 25  to  the  slotted  wall  results.  The 
results  are  now  seen  to  collapse  fairly  well  to  a common  curve  It  may  be  argued  that  if  the 
solid  wall  results  are  truly  interference  free,  then  such  a collapse  is  inevitable  from  the  method 
used  to  determine  the  value  of  / Nevertheless,  such  a result  is  a further  indication  of  the  validity 
of  the  corrections. 

In  summary,  the  application  of  the  complete  classical  linear  interference  theory  to  the 
slotted  wall  results  leads  to  fair  agreement  in  lift  curve  slopes,  indicating  a value  of  ti  /’  of  0-4. 

I inear  theory,  however,  badly  underestimates  the  magnitude  of  the  blockage  effects  leading  to 
corrected  Mach  numbers  in  error  by  up  to  0 -01.  The  introduction  of  a semi-empirical  blockage 
correction,  while  retaining  the  linear  theory  formulation  of  lift  interference,  leads  to  a similar 
scatter  in  the  corrected  lift  curve  slopes  for  a value  of  4 /’  of  0-51.  Fv  idence  from  several  sources 
indicates  that  such  a semi-empirical  approach  leads  to  corrected  data  which  are  very  nearly 
free  from  the  effects  of  wind  tunnel  wall  interference 


5.  COMPARISONS  WITH  OT1IF.R  DATA 
5.1  Comparisons  with  Other  Measurements 

A considerable  volume  of  test  data  has  become  available  in  recent  years  for  both  the 
aerofoil  sections  considered  in  the  present  tests  These  data  prov  ide  a useful  source  of  verification 
for  the  corrections  applied  to  the  data  in  the  present  tests  l or  the  BCik-l  section,  two  major 
sources  of  data  are  available  front  tests  conducted  at  the  NRC  in  Canada’  " and  from  the  ARA-’ 
in  the  United  kingdom  The  NRC  tests  were  conducted  in  the  NAI  15  60  two-dimensional 

test  facility,  with  the  top  and  bottom  walls  of  the  tunnel  having  6"„  porosity  for  the  earlier 
tests,  ’ and  20  - 5",,  porosity  for  the  later  tests  For  both  series  of  tests  the  model  chord  was  10' 
giving  a chord  to  height  ratio  (<  h)  of  0- 167.  No  corrections  for  blockage  (believed  to  be  small) 
or  lift  interference  have  been  applied  to  the  results.  I he  great  advantage  of  tilts  facility  is  its 
high  Reynolds  number  capability  (averaging  21  Id1’  for  the  present  tests),  but  this  is  obtained 
at  the  expense  of  aspect  ratio,  the  aspect  ratio  of  the  present  aerofoil  being  only  15,  Sidewall 
suction  is  applied  m the  vicinity  of  the  model  in  order  to  maintain  two-dimensional  flow  over 
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FIG.  25  COMPARISON  OF  PRESSURE  DISTRIBUTIONS  - BGK  V 

M = 0.705 
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FIG.  29  CORRECTED  SHOCK  POSITIONS 


(he  entire  model,  the  amount  of  applied  suction  being  determined  on  the  basis  of  flow  visual- 
ization and  experience  from  previous  tests  However,  there  remains  the  question  of  how  two- 
dimensional  the  flow  actually  is,  especially  in  the  presence  of  shock  waves  and  associated 
boundary-layer  separation. 

1 he  ARA  tests  were  conducted  in  its  IS'  S two-dimensional  facility  1 his  facility  has 
slotted  top  and  bottom  walls  giving  an  open  area  ratio  of  3 2"  Blockage  effects  in  this  tunnel 
are  considered  to  be  zero  but  lift  interference  corrections  were  applied  in  the  normal  manner 
with  003  and  Si  Oil.  The  model  tested  had  a chord  of  5"  and  a span  of  S'  giving 

a chord  to  height  ratio  of  0 278  and  an  aspect  ratio  of  1 6 The  Reynolds  number  based  on 
section  chord  was  maintained  at  6 10®  for  the  majority  of  the  tests,  dropping  to  4 10® 

at  a Mach  number  of  0 5 Unlike  the  NRC  tests  which  were  conducted  with  natural  transition, 
transition  was  fixed  on  both  surfaces  at  7 to  8"„  of  chord. 

figure  30  shows  a comparison  of  the  variation  of  lift  curve  slope  with  Mach  number  between 
the  corrected  data  or  the  present  tests,  and  the  results  from  the  NRC  and  ARA  tests  The  agree- 
ment on  the  whole  is  quite  good  with  the  ARA  results  falling  generally  above  those  of  the  present 
tests  as  would  be  expected  from  their  higher  test  Reynolds  number  The  NRC  results  would  be 
expected  to  lie  farther  above  again  on  the  basis  of  the  much  higher  Reynolds  number.  The 
lact  that  they  do  not  seems  to  indicate  either  a need  for  the  application  of  lift  interference 
corrections  in  a direction  applicable  to  an  open  jet,  or  a loss  of  lift  due  to  a lack  of  two- 
dimensionaliiy 

figure  31  shows  a comparison  ol  the  slotted  wall  design  pressure  distribution  for  the 
203  2 mm  chord  model  of  BGK-I  with  a similar  pressure  distribution  from  the  ARA  tests  This 
ARA  distribution  was  chosen  on  the  basis  that  it  exhibited  the  closest  agreement  with  the 
design  distribution.  The  Mach  numbers  (corrected  in  the  case  of  the  present  data)  agree  exactly, 
whereas  the  corrected  angles  of  incidence  and  lift  coefficients  are  somewhat  different  Whether 
these  differences  are  due  to  the  differences  in  test  Reynolds  number  or  to  incorrect  corrections 
for  the  effects  of  lift  interference  (in  one  or  both  cases)  is  at  this  stage  unclear.  Once  aga.n  a 
relative  loss  of  lift  in  the  upstream  region  of  the  lower  surface  is  apparent  for  the  present  203  2 n m 
tests,  indicating  the  presence  of  significant  streamline  curvature.  T his  effect  will  be  more  apparent 
in  this  case  since  the  magnitude  of  the  factor  5,  for  the  ARA  tests  is  similar  to  that  of  the  present 
tests  though  of  opposite  sign 

for  the  case  ol  the  NACA-0012  aerofoil  there  are  rather  more  recent  test  results  available 
for  comparison  Once  again  results  arc  available  from  the  N \1  tunnel  of  the  NRC.-®  For  these 
tests  a model  of  300  mm  chord  was  used,  giving  chord  Reynolds  numbers  up  to  40  x 10®. 
W ith  an  aspect  ratio  of  only  I - 27  the  effect  of  the  level  of  sidewall  suction  >n  the  two- 
dimensionality  ot  the  flow  was  quite  severe.  The  final  suction  rate  was  chosen  as  that  which  pro- 
duced a straight  separation  front  rather  than  using  the  normal  N AE  criterion  of  surface  stream- 
lines being  parallel  to  the  freestream  direction  over  the  complete  span.  The  justification  for  the 
use  ot  this  suction  rate  (about  half  the  rate  satisfying  their  normal  criterion)  is  stated  as  being 
"firstly  because  the  pressure  distribution  at  this  suction  level  appeared  to  match  pressure  data 
measured  at  ONI  RA.  and  secondly,  because  the  straight  as  opposed  to  the  curved  separation 
Iron!  was  telt  to  be  important  in  the  context  of  two-dimensional  simulation."-®  Tl.e  data  are 
considered  to  be  free  of  blockage  interference  but  corrections  have  been  applied  for  lift  interfer- 
ence based  on  both  theoretical  (classical  linear  theory  ) and  empirical  (streannvise  momentum 
balance)  grounds  The  results  quoted  here  are  limited  to  those  corrected  bv  the  application  of 
linear  theory. 

A very  useful  series  of  tests  were  conducted  on  the  NACA-0012  section  at  Calspan  in  1473.-' 
These  tests  were  conducted  in  the  Calspan  8'  8'  tunnel  with  slotted  top  and  bottom  walls 

giving  an  open  area  ratio  of  22  5",  The  model  had  a chord  of  6'  giving  a chord  to  height 
ratio  of  0 063.  and  the  results  would  therefore  be  expected  to  be  close  to  interference  free 
Calculations  with  linear  theory  at  Calspan  indicated  that  the  effects  of  blockage  would  amount 
to  less  than  0 I1’,,  and  those  of  lift  interference  to  less  than  2"„.  the  interference  being  of  the 
open  jet  type  The  model  had  a span  of  48  and  was  fitted  with  end  plates,  giving  a claimed 
effective  aspect  ratio  of  72  This  should  ensure  the  complete  absence  of  three-dimensional  effects 
No  corrections  were  applied  to  the  data,  anil  the  chord  Reynolds  number  was  maintained  at 
I x 10®. 

Further  results  are  available  from  tests  at  the  NI’l  -*  in  their  3b’  14'  transonic  tunnel. 
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FIG.  31  COMPARISON  OF  EXPERIMENTAL  DESIGN  PRESSURE  DISTRIBUTION:  - BGK 
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These  tests  were  conducted  with  a 10"  chord  model  (cjh  0-28)  of  14"  span  (aspect  ratio  I -4), 
The  tunnel  was  equipped  with  slotted  top  and  bottom  walls  with  an  open  area  ratio  of  33'',,. 
The  original  report28  claimed  that  the  conditions  for  these  tests  were  "close  to  those  giving 
blockage-free  and  lift-interference-free  results,  and  no  corrections  for  wall  constraint  have  been 
applied”.  Lately  however,  Collyer29  has  stated  that  this  is  now  believed  to  be  incorrect,  in  that 
although  there  is  substantial  experimental  evidence  to  suggest  that  the  conditions  are  indeed 
blockage-free,  lift  interference  is  significant,  and  gives  values  for  the  standard  lift-interference 
parameters  So  and  Si  ( 0 088  and  j 0 068  respectively).  The  test  Reynolds  number  varied  from 
I -7  x 10«  at  M 0-3  to  3 - 75  X I0«  at  M = 0 85. 

The  final  source  of  experimental  data  on  the  NACA-0012  is  some  French  tests  conducted 
by  ONF-RA  in  the  S3  tunnel  at  Modane  and  in  the  R1  tunnel  at  Chalais-Meudon.30  The  tests  in 
the  R1  tunnel,  with  working  section  dimensions  200  mm  70  mm,  were  conducted  with  porous 
top  and  bottom  walls  with  the  porosity  varying  between  0"o  (solid  walls)  and  l2-5°0.  Three 
models  were  tested  with  chords  of  60,  80  and  100  mm  respectively  (cjh  0-30,  0-40  and  0 50), 
with  all  models  spanning  the  70  mm  tunnel  dimension  (giving  aspect  ratios  of  1 17,  0-87  and 
0 -70).  From  comparisons  of  zero  lift  shock  positions  for  the  three  models  tested  in  walls  of  four 
different  porosities,  it  was  determined  that  the  walls  with  1 2 - 5 ‘ V,  porosity  gave  blockage  free 
conditions.  The  value  of  fi/P  calculated  from  linear  theory  to  be  relevant  to  the  condition  of 
zero  solid  blockage  (/3/P  1-28  for  perforated  walls)  was  taken  as  the  correct  value  for  the 

12 -5",,  porosity  walls  and  used  to  determine  lift-interference  parameters.  The  S3  tunnel  has 
test-section  dimensions  of  780  mm  x 560  mm  and  was  fitted  with  porous  top  and  bottom 
walls  of  9",,  porosity  for  these  tests.  The  model  tested  had  a chord  of  210  mm,  giving  a chord  to 
tunnel  height  ratio  of  0-27,  and  completely  spanned  the  tunnel  (giving  an  aspect  ratio  of  2 - 67). 
On  the  basis  of  equality  of  zero  incidence  drag  on  three  models  of  the  NACA-0012  section  of 
differing  chords,  it  had  been  determined  that  the  9%  porosity  gives  a condition  close  to  zero 
blockage.  Determination  of  porosity  factor  and  hence  lift-interference  parameters  had  not  been 
carried  out  at  the  time  of  the  tests.  The  Reynolds  number  for  the  tests  in  the  S3  tunnel  varied 
between  2 x I06  and  10  x 10s  and  in  the  Rl  tunnel  between  3 ■ 1 0,:  and  6 x I06. 

Figure  32  shows  a comparison  of  lift  curve  slopes  between  the  corrected  data  of  the  ARL 
tests  and  data  gathered  from  the  various  other  sources.  As  for  the  BGK-I  results,  the  overall 
agreement  is  quite  good,  but  several  points  merit  further  comment.  Firstly,  the  NPL  results, 
when  corrected  for  lift  interference  in  the  manner  suggested  by  Collyer,29  fall  well  away  from 
the  main  body  of  points,  whereas  the  uncorrected  data  are  in  general  agreement.  This  suggests 
that  the  statement  of  the  original  report28  that  the  results  were  close  to  being  free  of  lift 
interference  may  not  be  greatly  in  error.  Secondly  the  large  range  in  Reynolds  number  covered 
by  these  tests  (from  1 x 106  to  greater  than  40  x 1 0” ) does  not  produce  the  significant  differ- 
ences in  lift  curve  slope  that  might  be  expected.  Since  all  the  high  Reynolds  number  results  were 
obtained  from  tests  at  low  values  of  aspect  ratio,  loss  of  lift  due  to  three-dimensional  effects  is 
one  possible  explanation.  Finally,  it  is  noticeable  that  the  NRC  results  suffer  from  a distinct 
loss  of  lift  at  the  higher  Mach  numbers.  Again  this  can  probably  be  put  down  to  a lack  of  two- 
dimensionality  especially  in  the  region  aft  of  the  shock.  Indeed,  the  author  of  the  NRC  report,  in 
discussing  those  results,  suggests  that  it  may  perhaps  have  been  better  to  use  the  higher  suction 
rates,  corresponding  to  a more  two-dimensional  streamline  pattern  overall,  in  order  to  avoid 
this  problem.  He  comments  that,  in  comparison  with  the  ONF.RA  data,  the  NAE  measurements 
exhibit  a quicker  pressure  recovery  aft  of  the  shock,  leading  to  a distinct  loss  of  lift. 

Zero  lift  shock  positions  are  compared  in  Figure  33.  Again  the  agreement  between  the 
present  corrected  results  and  the  other  measurements  is  generally  quite  good.  The  expected  effect 
of  increasing  Reynolds  number  the  shock  location  moving  slowly  towards  the  trailing  edge 
is  not  generally  visible,  except  to  some  extent  in  the  NPL  results  at  higher  Mach  number  and  in 
the  opposited  direction  in  the  Calspan  results.  The  sensitivity26  of  shock  location  to  suction  rate 
in  the  NRC  tests  could  be  responsible  for  at  least  part  of  the  more  forward  shock  position,  but 
like  the  ONERA  Rl  tunnel  tests,  a general  lack  of  aspect  ratio  is  likely  to  be  the  major  factor. 
The  extremely  low  aspect  ratios  of  the  Rl  tunnel  tests  (I  ■ 17  to  0-70)  are  strongly  reflected  in  the 
extreme  forward  shock  locations. 
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FIG.  32  NACA  0012  LIFT  CURVE  SLOPE  - 

COMPARISON  WITH  OTHER  MEASUREMENTS 
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FIG.  33  NACA  0012  SHOCK  POSITION  - COMPARISON  WITH  OTHER  MEASUREMENTS 


5.2  Comparisons  with  Theory 


In  recent  years  there  have  been  impressive  advances  in  theoretical  methods  for  predicting 
the  pressure  distribution,  lift  and  drag  for  steady  inviscid  How  past  lilting  aerofoils  at  high  sub- 
sonic and  transonic  speeds.  This  development  began  with  the  proposal  by  Murman  and  Cole31 
of  a method  for  the  iterative  solution  of  the  transonic  small  perturbation  equations.  The  method 
was  further  developed  by  Garabedian  and  Korn4  3-  who  used  a conformal  transformation  due 
to  Sells33  combined  with  an  iterative  solution  of  the  exact  irrotational,  isentropie  equations  in  a 
manner  similar  to  that  used  by  Murman  and  Cole  It  was  recognized,  however,  that  in  the  range 
of  Mach  numbers  considered,  and  in  the  range  of  Reynolds  numbers  in  which  aerofoils  are 
designed  to  operate,  the  effects  of  viscosity  are  of  great  importance,  and  hence  several  attempts 
have  recently  been  made  to  develop  theoretical  models  which  include  these  effects. -s-31  30  Of  these 
methods,  only  those  developed  at  the  RAF-9-38  make  any  attempt  to  correctly  take  account  of  the 
effects  of  the  thickness  and  curvature  of  the  aerofoil  wake,  as  well  as  the  displacement  effect 
of  the  aerofoil  boundary  layers. 

The  verification  of  these  various  theoretical  methods  still  awaits  the  provision  of  experi- 
mental measurements  which  are  truly  two-dimensional  and  free  from  the  effects  of  wind-tunnel 
wall  interference,  and  this  was  indeed  one  of  the  reasons  for  the  initiation  of  the  CAARC 
program.3  Computer  programs  are  now  readily  available  lor  calculating  inviscid  flows,  but 
unfortunately  this  is  not  yet  the  case  for  the  more  recently  developed  viscous  theories.  Com- 
parisons of  these  theories  with  experimental  results  are  therefore  limited  to  the  relatively  few 
cases  which  have  been  reported  in  the  literature. 

Figure  34  compares  corrected  lift  curve  slope  data  from  the  present  tests  with  inviscid 
theory.  As  expected,  the  inviscid  theory  greatly  over-estimates  lilt  curve  slope,  being  approxi- 
mately 20",,  greater  than  the  corrected  data  throughout  the  range  of  Mach  number  Unfortunately 
no  comparable  results  are  available  from  viscous  theories  for  this  section.  It  ae  case  of  the 
NACA-0012  section,  however,  results  of  theoretical  viscous  calculations  tire  available  from  the 
methods  of  Firmin36  (quoted  by  Lock3)  and  Collyer,--'  and  these,  together  with  an  inviscid 
solution,  are  compared  with  the  corrected  data  in  f igure  35.  Once  again,  the  inviscid  solution 
overestimates  the  lift  curve  slope,  in  this  case  by  an  average  of  25",,.  the  difference  increasing 
slightly  with  Mach  number.  With  the  effects  of  viscosity  included,  the  difference  between  theory 
and  experiment  is  reduced  to  about  15",,.  The  results  from  (.'Oliver's  theory,  which  would  be 
expected  to  be  more  accurate  (due  to  its  utilization  of  the  exact  equations  to  represent  the 
inviscid  part  of  the  flow)  than  that  of  firmin  (which  uses  the  transonic  small  perturbation 
equations),  show  the  greater  difference.  However,  the  results  shown  in  I igurc  35  are  for  a 
Reynolds  number  of  3-5  • 10®  compared  to  an  average  Rey  nolds  number  of  the  experimental 
results  of  I 2 x I0fi.  The  difference  in  lift  curve  slope  due  to  this  change  in  Reynolds  number, 
on  the  basis  of  published  results3  from  I irmin's  theory,  would  amount  to  about  4 In  fact, 
correcting  Collyer's  results  in  this  manner  produces  quite  good  agreement  between  the  two 
theories,  but  there  still  remains  a difference  of  more  than  10",,  between  theory  and  the  corrected 
experimental  results. 

Figure  36  compares  zero  lift  shock  positions  predicted  by  inviscid  theory  for  the  NACA- 
0012  with  those  of  the  corrected  experimental  results.  Also  included  is  a single  point  available 
from  Collyer's  theory.  The  agreement  between  theory  and  experiment  is  again  quite  good 
This  type  of  comparison  is,  however,  not  strictly  valid.  Hie  governing  equations  used  in  the 
inviscid  theory  (and  in  the  viscous  theory  for  the  inviscid  part  of  the  llow)  arc  isentropie  and 
irrotational,  and  hence  the  compression  jumps  which  evolve  as  approximations  to  shock  waves, 
are  only  approximations  to  the  true  Rankine-Hugoniot  shock  jumps,  1 urthermore.  shock  waves 
are  smeared  over  several  mesh  spacings.  and  differences  in  shock  positions  between  non- 
conservative schemes  (as  used  here)  and  the  fully  conservative  relaxation  schemes  are  significant 
(See  for  example  Bailey37  or  Murman38  for  their  discussion  of  these  points  ) 

Among  the  published  results  for  the  BGk-l  section  from  Collyer's  theory,  there  arc  for- 
tunately several  cases  in  which  Mach  numbers  and  angles  of  incidence  are  sufficiently  close  to 
the  corrected  experimental  values  for  direct  comparison  of  pressure  distributions  fhree  cases 
are  shown  in  Figures  37  to  39.  The  first  is  the  experimentally  determined  .lesign  point  for  the 
203-2  mm  model.  The  second  is  for  the  same  Mach  number,  but  at  higher  incidence,  and  shows 
the  appearance  of  a fairly  strong  shock,  finally,  at  an  incidence  lower  than  the  design  value. 
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FIG.  34  BGK-1  LIFT  CURVE  SLOPE  - COMPARISON  WITH  THEORY 


Inviscid  theory  (Garabedian  — Ref.  31) 


IG.  35  NACA  0012  LIFT  CURVE  SLOPE  - COMPARISON  WITH  THEORY 


FIG.  36  NACA  0012  SHOCK  POSITION  - COMPARISON  WITH  THEORY 
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INVISCID  THEORY  (BAUER  ET.AL.  - REF.  32) 

M=  .760  AL=  0.50  CL=  0.757  CM=-0.171  CD=  0.011 

VISCOUS  THEORY  (COLLYER  - REF. 29) 

M=  .760  AL=  0.50  C!M=  0.554  R=  6.000 

+ X EXPERIMENT  203. ZMH  CHORD  SLOTTED  WALLS  CORRECTED 
M=  -760  AL=  0.50  CN=  0.533  CM=-0.1Z9  R=  1.646 


FIG.  37  COMPARISONS  OF  EXPERIMENTAL  AND  THEORETICAL  PRESSURE 
DISTRIBUTIONS  - BGK-1;  M = 0.76  a = 0.5 
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INVISCID  THEORY  (BAUER  ET  .AL . - REF.3Z) 


\t=  .758  AL=  1 .60  CL=  0.997  CM=-0.204  CDt  0.033 

VISCOUS  THEORY  (COLLYER  - REF.  79) 

M=:  .760  AL=  1 .60  CN=  0.744  R=  6.000 

+ X EXPERIMENT  101 .6MM  CHORD  SLOTTED  WALLS  CORRECTED 
M=  .758  AL-  1 .59  CN=  0.69Z  CM=-0.1ZZ  R=  0.794 


FIG.  38  COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL  PRESSURE 
DISTRIBUTIONS  - BGK  1,  M = 0.76;  a = 1.5° 
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INVISCID  THEORY  (BAUER  ET  .AL . - REF.  32) 

M=  -756  AL=  0.(30  CL=  0.630  CM=-0.155  CD=  0.004 


VISCOUS  THEORY  (COLLYER  - REF.  29) 

M=  .760  AL=  0.00  CN=  0.444 


R=  6.000 


+ X EXPERIMENT  101 .6M^  CHORD  SLOTTED  WALLS  CORRECTED 
M=  .758  AL=— 0 .06  CN=  0.397  CM=-0.1ZS  R=  0.794 


FIG.  39  COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL  PRESSURE 
DISTRIBUTIONS  - BGK  1;  M = 0.76,  a = 0° 


but  Still  at  the  same  Mach  number  we  see  the  typical  double  shock  system  predicted  by  both 
inviscid  and  viscous  theory  and  found  to  some  extent  in  experiment.  Also  included  on  these 
figures  are  inviscid  calculations  for  the  same  Mach  numbers  and  angles  of  incidence.  In  all 
cases  the  agreement  between  the  viscous  theory  and  experiment  in  better  than  that  of  the  inviscid 
theory.  Considering  the  difference  in  Reynolds  number  between  viscous  theory  and  experiment, 
the  overall  agreement  is  excellent,  the  theory  predicting  accurately  all  the  major  feature  of  the 
pressure  distributions.  However,  lift  coefficients  are  still  overestimated  by  up  to  10",,,  the  closest 
agreement  occurring  somewhat  surprisingly  at  the  experimental  design  condtiion,  where  the 
difference  is  less  than  50o. 

6.  DISCUSSION 

From  the  comparisons  of  the  previous  sections,  it  is  evident  that  the  application  of  correc- 
tions based  on  classical  linear  interference  theory,  with  the  addition  of  a semi-empirical  formula- 
tion of  the  effects  of  blockage  in  the  slotted-wull  case,  produces  data  which  form  a self-consistent 
set,  as  well  as  being  in  general  agreement  with  data  from  other  sources.  Several  significant 
points  stand  out.  Firstly,  the  corrections  continue  to  produce  good  agreement  ever,  when  the 
Mach  number  is  such  as  to  produce  significantly  supercritical  flow  over  the  aerofoil  upper 
surface.  Since  the  basis  from  which  linear  theory  is  derived  is  certainly  not  valid  under  super- 
critical conditions,  this  result  is  somewhat  surprising  The  present  tests  do  indicate  however, 
that  at  least  up  to  Mach  numbers  at  which  lift  coefficients  begin  to  decrease  due  to  shock 
induced  separation  on  the  upper  surface  and  the  appearance  of  a supersonic  /one  on  the  lower 
surface,  the  application  of  linear  theory  still  produces  quite  good  results  Secondly,  the  basic 
linearized  approach  to  the  problem  inherent  in  classical  theory,  which  allow  s the  modelling  of  the 
situation  by  three  independent  corrections  for  solid  and  wake  blockage  and  lift  interference,  is 
generally  supported.  The  overall  implications  of  this  approach  are  that  if  corrections  are  applied 
independently  to  the  Mach  number  and  angle  of  incidence,  then  (he  pressure  distributions 
obtained  for  identical  corrected  conditions  should  also  agree  (If  pressure  distributions  are 
compared  in  the  form  of  pressure  coefficients  rather  than  p.essure  ratios,  then  care  must  be  taken 
to  correct  the  pressure  coefficients  for  the  effects  of  blockage  on  dynamic  pressure.  It  is  for  this 
reason  that  all  pressure  distributions  considered  in  this  report  are  expressed  in  the  form  of 
pressure  ratios.)  This  has  indeed  been  found  to  be  the  case  in  comparisons  between  sets  of 
corrected  data  from  the  present  tests  and  also  in  comparisons  with  corrected  data  from  other 
sources.  This  type  of  comparison  is  of  course  only  valid  so  long  as  the  effects  of  streamline 
curvature  remain  small.  As  noted  earlier,  the  presence  of  streamline  curvature  leads  to  the 
implication  that  a warped  model  should  be  tested  to  simulate  the  uniform  flow  situation,  the 
degree  of  warpage  being  a function  of  angle  of  incidence.  In  the  present  tests,  the  maximum 
correction  to  lift  coefficient  due  to  the  effects  of  streamline  curvature  amounts  to  less  than  3'\, 
in  the  solid  wall  case  and  about  I for  the  slotted  walls  Thus  the  streamline  curvature  effects 
in  the  present  tests  remain  small  enough  not  to  seriously  affect  the  validity  of  the  corrections, 
although  some  differences  due  to  such  effects  are  visible  in  pressure  distributions  of  the  larger 
(203  2 mm  chord)  models  of  both  sections  Finally,  it  is  worth  noting  that  the  scatter  in  the 
agreement  of  the  corrected  solid  wall  results,  and  in  the  comparison  of  corrected  results  from 
both  solid  and  slotted  wall  tests  is  less  than  might  be  expected.  As  stated  previously  the  errors 
in  measured  lift  curve  slopes  are  expected  to  be  on  the  order  of  2 , ...The  fact  that  scatter 
between  corrected  results  is  only  slightly  greater  than  this  level  is  indeed  quite  remarkable 

Despite  the  success  of  classical  linear  theory  in  correcting  the  present  results,  its  application 
is  not  without  problems  The  major  problem  area  is  of  course  the  inability  of  classical  theory  to 
correctly  represent  the  effect  of  blockage  in  the  slotted  wall  situation.  It'  predictions  have  been 
seen  to  underestimate  the  required  corrections  to  Mach  number  by  up  to  a factor  of  live  Several 
possible  mechanisms  to  explain  this  anomaly  have  been  investigated  Frrors  in  test  section  total 
pressure  measurement,  which  for  the  present  tests  was  determined  from  the  static  pressure  in  the 
entrance  to  the  contraction,  are  possible  due  to  the  finite  velocity  in  the  contraction  However, 
the  large  contraction  ratio  of  the  \RI  transonic  tunnel  (22  I)  minimizes  this  source  of  error, 
the  error  amounting  to  less  than  0 0003  in  Mach  number  at  a Mach  number  of  O S.  A further 
possibility  is  the  suggestion  by  Read  ri  al  M that  the  introduction  of  models  with  large  (greater 
than  2'’,,)  blockage  ratios  could  reduce  plenum  chamber  pressure  below  the  calibrated,  empty- 
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tunnel  values  at  high  subsonic  Mach  numbers  Such  an  action  would  have  the  effect  of  in- 
creasing the  Mach  number  calculated  Iront  static  pres  ure  in  the  plenum  chamber  above  that 
without  the  model  present  Such  effects  were  found  b>  Parker'  when  studying  the  flow  generation 
properties  of  a set  of  slotted  walls  in  the  Propulsion  Wind  I unnel  I aeility  at  the  Arnold  Fngin- 
eertng  Development  C entre  The  magnitude  ot  this  effect,  from  Parker's  measurements,  appeared 
to  increase  with  Mach  number,  reaching  almost  0 01  in  Mach  number  at  a Mach  number  of  0 9. 
A short  investigation  was  undertaken  in  the  API  transonic  tunnel  to  discover  whether  such  an 
effect  was  the  cause  of  the  discrepencv  in  blockage.  Centreline  side-wall  pressures  were 
measured  upstream  of  the  model  position  for  the  empty  tunnel  condition,  and  with  the 
203  2 mm  chord  N At  A-0012  model  installed  in  the  tunnel  at  an  angle  of  incidence  of  3 . No 
difference  in  the  relationship  between  plenum  chamber  static  pressure  and  upstream  test- 
section  wall  pressures  was  evident  between  the  empty  tunnel  and  model  installed  cases.  During 
this  investigation  a check  was  also  made  of  the  possibility  of  the  plenum  chamber  static  pressure 
varying  with  position  in  the  plenum  chamber  due  to  highly  curved  (lows  through  the  slotted 
walls.  Once  again,  however,  no  detectable  differences  were  apparent 

It  has  been  suggested  (e  g.  see  Lock3)  that  unlike  the  predictions  of  linear  theory,  the  ratio 
P and  hence  8U  and  &i  will  probably  be  a function  of  Mach  number.  The  slightly  increased 
amount  of  scatter  towards  the  high  Mach  number  range  in  the  corrected  lift  curve  slope  results 
of  the  present  tests  would  tend  to  support  this  proposal.  However,  in  this  region,  where  the  (low 
is  supercritical,  linear  theory  no  longer  applies,  and  allowing  ; P to  vary  with  Mach  number 
seems  to  be  simply  an  artificial  attempt  to  extend  the  range  of  validity  of  linear  theory  . 

As  well  as  being  subject  to  the  effects  ol  wind  tunnel  wall  interference,  transonic  wind  tunnel 
tests  on  two-dimensional  aerofoils  are  also  very  sensitive  to  the  magnitude  of  test  Reynolds 
numbers.  The  present  tests  provide  a good  opportunity  for  studying  Reynolds  number  effects 
both  in  the  tests  themselves  and  in  comparisons  with  other  data,  both  experimental  and 
theoretical.  Tests  in  the  slotted  walls  on  the  203  2 mm  chord  models  of  both  sections  at  the 
Reynolds  number  of  the  ltd  -6  mm  chord  models  indicated  that  the  effect  of  this  halving  of 
Reynolds  number  on  lift  curve  slopes  was  quite  small  of  the  same  order  as  the  scatter  (about 
d 2 “ ) in  the  data.  The  lower  Reynolds  number  data,  however,  did  show  a noticeably  lower 
trend,  as  would  be  expected.  In  the  corrected  data,  for  both  solid  and  slotted  walls,  the  101  -6 
mm  chord  results,  although  within  the  general  scatter,  do  produce  noticeably  lower  value's  of 
lift  curve  slope.  This  is  most  noticeably  in  the  N At  A-0012  results  from  the  slotted  walls  shown 
in  Figure  33.  Whether  this  effect  is  due  to  the  difference  in  Reynolds  numbers,  or  to  an  inability 
of  linear  theory  to  correctly  account  for  differences  in  model  si/e.  cannot  be  determined.  It  is 
however,  interesting  to  note  that  the  results  of  I it  min's  "'  theory  predict  a change  in  lift  curve 
slope  ot  just  under  I -5°,,  for  this  difference  in  Rey  nolds  n nber  As  noted  earlier,  the  expected 
differences  in  lift  curve  slope  due  to  significant  variations  in  Reynolds  number  were  not  found 
in  comparisons  of  the  present  data  with  other  measurements  Again,  based  on  Firmin's  theory, 
the  difference  in  lift  curve  slope  between  the  lowest  (I  10s)  and  highest  (40  10s)  Reynolds 

number  data  would  amount  to  1 1"„  at  a Mach  number  of  0-7.  Since  all  the  available  data  at 
Reynolds  numbers  greater  than  those  of  the  present  tests  were  obtained  on  models  with  aspect 
ratios  less  than  I -6.  the  appearance  of  three-dimensional  effects  would  seem  to  be  a likelv 
cause  of  the  general  lack  of  lilt  It  would  appear  that  in  their  quest  for  the  highest  possible 
Reynolds  number,  model  designers  have  erreu  towards  aspect  ratios  which  are  too  small  to 
produce  truly  two-dimensional  data. 

Comparisons  of  viscous  theory  w ith  the  corrected  results  indicate  that  these  theories  have 
reached  a stage  where  they  can  predict  the  overall  features  of  the  flow  over  aerofoils  at  transonic 
speeds  to  a high  degree  of  accuracy.  There  remains,  however,  a general  tendency  for  all  theories 
to  overestimate  lift  coefficients,  the  error  amounting  to  some  10  . for  the  best  of  the  theories. 
This  is  even  true  of  the  method  due  to  Collyer  which  has  the  most  complete  representation  of  the 
physics  of  the  situation,  including  contributions  from  wake  displacement  and  curvature  effects. 

Betore  completing  this  discussion  of  the  application  of  classical  interference  theory,  the  effect 
on  the  quality  of  the  data  of  testing  in  wind  tunnels  equipped  with  walls  of  different  open  area 
ratios  should  be  considered.  As  shown  previously . the  slotted  w alls  used  in  the  present  tests,  whose 
open  area  ratio  of  16-5",,  was  selected  primarily  on  the  basis  of  achieving  supersonic  operation 
through  diffuser  suction,  require  large  corrections  for  blockage,  amounting  to  about  0 01 5 
in  Mach  number  at  ,\1  0-8.  To  obtain  blockage  free  conditions,  it  has  been  suggested  that 
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milk'll  smaller  open  arc.i  ratios  arc  required.  in  I he  range  to  -I  liu  slotted  walls  However.  the 
open  walls  ol  the  present  tests  ih>  have  the  advantage  of  requiring  small  streamline  eui\atu;e 
eorreetions  As  the  open  area  ratio  is  decreased  the  \ahie  ol  Ai  ( 0 Os  t,,i  the  present  walls) 

will  ileerease  to  zero,  a ml  then  increase  again  (with  opposite  sign)  tow.uds  the  closed  wall  value 
of  it  ’(  0 1.1)  funnels  in  which  the  walls  have  been  optimized  towauls  blockage  Tree  conditions 

(lot  example  the  IS  S'  \KA  tunnel,  with  an  open  men  i.itio  ol  1 1 „)  appear  to  sutler 
from  a larger  value  ol  <S|.  tending  towauls  the  closed  wall  value  ( O il  Cot  the  AKA  tunnel) 
ll  the  eflects  ol  blockage  can  be  successfully  calculated  liom  lineal  theotv,  ot  computed  b\ 
empirical  means  within  the  overall  framework  ol  linear  theory  as  has  been  done  here,  the 
absence  of  blockage  free  conditions  would  seem  to  be  a small  price  to  pay  lor  the  mininii/ation 
ol  uncorreclable  streamline  curvatuie  eflects  faking  tins  appioaclt  to  its  logual  conchisnvn 
would  k»l  course  result  in  the  use  ol  walls  with  open  men  talios  aiiauged  such  that  streamline 
curvature  etVects  did  not  exist  (fv,  0)  in  genei.il.  such  walls  would  not  produce  blockage  free 
conditions,  but  this  type  ol  approach  mav  lead  to  bet  lei  overall  accuiacv  for  two-dimensional 
aerofoil  testing  at  high  subsonic  speeds 

7.  t ()M  I l SIONS 

I he  results  ol  transonic  wind  tunnel  tests  on  two  khniensioii.il  aerofoils  have  been  analvscvl 
w illv  respect  to  the  etVects  of  wind  tunnel  wall  interference  fltese  tests  vveie  comlucted  on  two 
geometrical!!  similar  models  >>1  each  of  two  aeioloil  sections  the  N At  \ POI.'  and  B(iK  I 
sections  and  covered  a range  in  Mach  number  liom  0 s tv'  0 S.'  Model  chords  dilVered  by  a 
lactoi  of  (wo.  giving  model  chord  to  tunnel  height  ratios  of  0 ,'s  and  (i  I.’'  In  order  to 
examine  the  validity  «'l  classical  subsonic  iinean/Ckl  wall  inlet leience  theoiv  without  the  com 
plications  id  the  doubtful  representation  of  non  solid  wall  bouiulmv  conditions,  the  results  of 
tests  conducted  in  test  sections  lilted  with  both  solid  and  slotted  walls  were  considered  I he 
results  have  been  corrected  using  the  predictions  of  classical  linear  theoiv  I oi  the  solid  wall 
tests,  these  corrections  appeal  to  produce  icstills  which  aie  veiv  dose  to  being  free  of  the  effects 
ot  wall  interference  In  the  ease  ot  the  slotted  wall  tests,  hovvcvei.  I r ne.it  theoiv  significantly 
underestimated  the  degree  of  blockage  evident  in  the  data  I 'sing  the  foim  ol  the  variation  of 
blockage  with  Mach  number  suggested  bv  lineal  theory,  and  an  empirical  constant  devised  from 
the  test  data,  a satisfactory  representation  of  the  blockage  collection  was  obtained  fhis 
empirical  blockage  correction  was  then  combined  with  the  variation  of  lift  interference  predicted 
Ironi  line  n theory . ami  the  crossflow  characteristics  (,•’  /’)  ,>|  ihe  slottevl  walls  were  then  obtained 
bv  comparison  of  the  slotted  wall  data  with  the  collected  solid  wall  data  (taken  as  being  free  of 
interference)  I his  procedure  led  to  a value  oft)  x|  tot  /•  and  values  of  lift  inteifei-nce  para 
meters  .\i  ami  tX ( » *1  0 I b ami  0 ()>  res|veclivelv 

Comparisons  between  the  sets  ot  corrected  data,  and  with  data  from  other  sources, 
indicated  that  the  correction  scheme  thus  derived  did  indeed  produce  results  which  were  largclv 
lice  ol  wall  inlet lerence  elfects  Such  cotnpat tsons  als<>  venlieil  that  the  simple  uulepeudeut 
corrections  1 1 » Mach  nunibei  ami  angle  «*l  incidence  implied  bv  lineal  theory,  ih>  ptoduce  results 
which  lepresent  the  same  llow  conditions  I lie  appeal  .nice  ol  siieanihne  curvature  elfects,  which 
invalidate  such  independent  collections,  was  limited  to  the  largei  choid  (.'(H  .'mm)  nii'ilels. 
ami  at  all  tunes  lem.uneil  small  Ihe  implications  I'm  future  two-dimensional  testing  in  the 
slottevl  wall  test  section  arc  therefore  as  follows  1 o avoid  am  streamline  curvature  elfects. 
model  chords  should  be  somewhat  stttallci  than  those  ol  the  largei  models  tested  here  I o 
retain  ease  ot  mamif.ictme  and  to  allow  installation  of  suthcicnt  piessure  ontices  to  adequate)! 
deline  the  sm lace  pressure  ihstribution.  a chord  i>l  about  150mm  wimlil  appeal  oplinnim 
I lie  application  i'l  the  semi-empirical  collection  scheme  ileriveil  lieie  shouhl  then  produce 
results  which  are  free  ol  wind  tunnel  wall  interference 

In  eontparing  the  present  corrected  lesiilts  with  data  liom  othei  sv*utces,  several  interesting 
features  were  apparent  I itstlv.  among  most  modern  tests,  especially  those  conducted  at  highei 
Reynolds  numbers,  there  appears  to  be  an  unloitiinate  tendenev  to  saciitice  aspect  ratio  in  the 
search  foi  the  highest  possible  Reynolds  uumhet  I his  has  been  taken  to  such  an  extent  that 
aspect  ratios  ol  little  more  than  unit \ have  become  the  norm  I loin  conipni isons  i>l  the  present 
data  with  othei  experimental  measurements  it  was  evident  that  mam  of  the  higliet  Remolds 
nunibei  tests  ihd  not  show  the  expected  highei  values  »»t  hit  cmve  slo|ve,  and  it  is  likely  that  the 
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lack  ol  aspect  ratio  in  these  tests  is  the  prime  cause.  Secondly,  comparisons  with  the  latest 
theoretical  methods  for  calculating  the  viscous  flow  around  an  aerofoil  at  transonic  speeds 
indicated  that  such  theories  are  now  capable  of  fairly  accurate  predictions.  Major  features  of 
pressure  distributions  are  represented  quite  well,  even  in  areas  expected  to  be  sensitive  to  small 

changes  in  flow  conditions.  There  remains,  however,  a general  tendency  to  overestimate  lift 
coefficients  by  about  10°,,. 


i 

* 


I 


REFERENCES 

1.  Report  of  Meeting  No.  1!  of  CAARC,  held  in  Australia.  17  26  September  1473  C C 721 

2.  Moss.  G.  F.  Report  on  meeting  of  Co-ordinators  in  the  Field  of  Aerodynamics,  held  in 
India,  24  November  to  4 December  1472.  C.C.  704.  January  1474. 

.V  Lock.  R C.  A proposal  for  co-operative  aerofoil  tests  in  Commonwealth  transonic  wind 
tunnels.  C.C.  724.  March  1475. 

4 Bauer,  F.,  Garabcdian,  P.,  and  Korn.  D.  G.  Supercritical  wing  sections.  (Lecture  notes  in 
economics  and  mathematical  systems,  Vol  66  ) Spiingcr-Verlag.  Berlin,  1472. 

5.  kaeprzynski,  J.  J.,  Ohman,  L IT.  Garabcdian,  P.  R . and  Korn.  D.  G Analysis  of  the  flow 
past  a shockless  lifting  aerofoil  in  design  and  otT-dcsign  conditions  NRC.  NAF  Report 
LR-554,  1471. 

6.  Kaeprzynski,  J.  J.  An  experimental  analysis  and  buffet  investigation  of  the  shockless  lifting 
aerofoil  No.  I.  NRC.  NAF  Report  LR-569.  1473. 

7.  Pollock,  N.,  and  Fairlie,  B D.  An  investigation  of  supercritical  aerofoil  BGK-I  Part  I. 
Near  design  point  tests  and  comparisons  with  theory  ARL  Report  A 144.  1475. 

8.  Pollock.  N.,  and  Fairlie,  B.  D.  An  investigation  of  supercritical  aerofoil  BGK-I  Part  2: 
Test  Reynolds  number  requirements  and  transition  fixing.  ARL  Report  A 145.  1475. 

4.  Fairlie.  B.  D.,  and  Pollock.  N.  An  investigation  of  supercritical  aerofoil  BGK-I  Part  3: 
Extended  subsonic  and  transonic  tests  ARL  Report  A 146.  1475. 

10.  Pollock.  N.,  and  Fairlie.  B.  D.  Transonic  wind  tunnel  tests  on  a NACA  (XII 2 aerofoil.  ARL 
Report  A/ 148.  1476. 

1 1 Garner.  11  C„  Rogers,  E.  W.  E..  Acum,  W.  1 A , and  Maskcll.  F.  C.  Subsonic  wind-tunnel 
wall  corrections.  AGARDograph  104,  1466. 

12.  Berndt.  S.  B.,  and  Sorensen,  H.  Flow  properties  of  slotted  wails  for  transonic  test  sections. 
AGARD  CP-174.  1476. 

13.  Woods,  L.  C.  On  the  theory  of  two-dimensional  wind  tunnels  with  porous  walls.  Proc.  Roy 
Soc.  A.  233.  74.  1455. 

14  Fairlie,  B.  D , and  Pollock.  N.  Transonic  wind  tunnel  tests  on  a series  of  two-dimensional 
aerofoils  in  a solid  wall  test  section.  ARL  Aero.  Note  A 384.  1474. 

15.  Baldwin.  B.  S..  Turner.  J.  If.  and  Knechtal.  E.  I).  Wall  interference  in  wind  tunnels  with 
slotted  and  porous  boundaries  at  subsonic  speeds.  NACA  TN  3176.  1454 

16.  Wright.  R H.  The  effectiveness  of  the  transonic  wind-tunnel  as  a device  for  minimising 
tunnel-boundary  interference  for  model  tests  at  transonic  speeds.  AGARD  Rcpt.  244.  1454 

17.  Catherall.  D.  On  the  evaluation  of  wall  interference  in  two-dimensional  ventilated  wind 
tunnels  by  subsonic  linear  theory  RAF  TR76I34  1476. 

18.  Pindzola,  M , and  Lo.  C.  F Boundary  interference  at  subsonic  speeds  in  wind  tunnels  with 
ventilated  walls  AEDC-TR-64-47.  1464. 

14.  Davis,  D.  D . and  Moore.  D.  Analytical  study  of  blockage-  and  lift-interference  corrections 
for  slotted  tunnels  obtained  by  the  substitution  of  an  equivalent  homogeneous  boundarv 
for  the  discrete  slots.  NACA  RM  L53E07B.  1453. 


r 


i 


< 


20.  Barnwell,  R W.  Improvements  in  the  slotted-wall  boundary  condition.  Proc.  AIAA  9th 
Aero.  Testing  Conf.,  Texas,  1976. 

21.  Chen,  C.  F.,  and  Mears.  J.  W.  Experimental  and  theoretical  study  of  mean  boundary  con- 
ditions at  perforated  and  longitudinally  slotted  wind  tunnel  walls  AIT)C  TR-57-20,  1957. 

22.  Lo,  C.  H.,  and  Oliver,  R.  H.  Boundary  interference  in  a rectangular  wind  tunnel  with 
perforated  walls.  AEDC  TR-70-67,  1970. 

23.  Monti.  R.  Wall  corrections  for  airplanes  with  lift  in  transonic  wind  tunnel  tests.  AGARD 
Advisory  Report  No.  36,  1971. 

24.  Pearcey,  H.  H.,  Sinnott,  C.  S.,  and  Osborne,  J.  Some  effects  of  wind  tunnel  interference 
observed  in  tests  on  two-dimensional  aerofoils  at  high  subsonic  and  transonic  speeds. 
AGARD  Report  296.  1959 

25.  Hammond,  B.  F.  L.  Results  of  tests  on  aerofoil  M.8I  5 (BGK-I)  in  the  ARA  two-dimensional 
tunnel.  ARA  Model  Test  Note  M.81/5,  1978. 

25.  Ramaswamy,  M.  A.,  Krishnamurthy,  K , and  Ramachandra  Sharma.  Transonic  wind  tunnel 
tests  on  a supercritical  aerofoil.  NAI.  TM  AE-TM-24-74,  1974. 

26.  Dixon,  R.  C.  High  Reynolds  number  investigation  of  an  ONI  RA  model  of  the  NACA  0012 
airfoil  section.  NRC  LTR-H A-5X5/0069,  1975. 

27.  Vidal.  R.  J.,  Caltin,  P.  A.,  and  Chudyk.  D.  VV.  Two-dimensional  subsonic  experiments  with 
an  NACA-0012  airfoil.  Calspan  Rept.  No.  RK-5070-A-3,  1973. 

28.  Gregory,  N.,  and  Wilby,  P.  G.  NPL96I5  and  NACA-0012.  A comparison  of  aerodynamic 
data.  ARC  CP  No.  1261,  1973. 

29.  Collyer,  M.  R.  An  extension  to  the  method  of  Garabcdian  and  Korn  for  the  calculation  of 
transonic  flow  past  an  aerofoil  to  include  the  effects  of  a boundary  layer  and  wake.  RAF 
TR77I04,  1977. 

30.  Bazin.  M.,  Bernard-Guelle,  R.,  and  Pontezierc,  J.  Critique  des  techniques  d'essais  de  ptoltls 
transsoniques.  Premiere  partie:  Dispositif  d'essais  industriels  a S3MA.  L'Aeronautiquc  et 
L'Astronautique,  31,  1971. 

Deuxiemc  partie:  fitude  experimentale  des  corrections  de  parois  a Rl  Ch.  L'Adro- 
nautique  et  L'Astronautique,  32,  1971. 

31.  Murman,  E.  M.,  and  Cole.  .1.  D.  Calculation  of  plane  steady  transonic  flow.  AIAA  J,  9. 

1 14  21.  1971. 

32.  Bauer,  F.,  Garabcdian,  P.  R.,  Korn,  D.  G.,  and  Jameson  A.  Supercritical  wing  sections  II 
Lecture  Notes  in  Economics  and  Mathematical  Systems.  Springer- Verlag.  1975. 

33.  Sells,  C.  C.  L.  Plane  subcritical  flow  past  a lifting  aerofoil.  Proc.  Roy.  Soc.  A.  308,  377  401. 
1968. 

34.  Bavitz,  P.  C.  An  analysis  method  for  two-dimensional  transonic  vicsous  flow  NASA  TN 
D-7718,  1975. 

35.  Bauer,  F.,  and  Garabcdian.  P.  R.  Computer  simulations  of  shock  wave  boundary -layer 
interactions.  Comm.  Pure  and  App.  Maths.  36.  659  65.  1973. 

36.  Firmin,  M.  C.  P.  Calculation  of  the  pressure  distribution,  lift  and  drag  on  aerofoils  at  sub- 
critical  conditions.  Part  I Interim  method.  RAF  I R 72235.  1973. 

37.  Bailey,  F R On  the  computation  of  two-  and  three-dimensional  steady  transonic  flows  by 
relaxation  methods.  Prog,  in  Num.  Fluid  Dyn.,  von  Kantian  Inst,  for  Fluid  Dyn„  Lecture 
Series  63,  1974. 

38.  Murman,  F.  M.  Analysis  of  embedded  shock  waves  calculated  by  relaxation  methods 
Proc.  AIAA  Comp.  Fluid  Dyn.  Conf.,  Palm  Springs,  1973. 


L 

* 


I 


39  Recd-  T D • 1 C • and  Cooksey.  J.  M.  Calibration  of  transonic  and  supersonic  wind 

tunnels.  NASA-CR-2920.  1977. 

40.  Parker.  R.  L.  Flow  generation  properties  of  live  transonic  wind  tunnel  test  section  wall 
configurations.  AEDC-TR-75-73,  1975. 


DOCUMENT  CONTROL  DATA  SHEET 


Security  classification  of  this  page:  Unclassified 


1.  Document  Numbers 

2.  Security  classification 

(a)  AR  Number: 

(a)  Complete  document  : 

AR  -001-596 

Unclassified 

(b)  Document  Series  and  Number: 

(b)  Title  in  isolation: 

Aerodynamics  Report  151 

Unclassified 

(c)  Report  Number: 

(c)  Summary  in  isolation: 

ARL — Aero-Report  151 

Unclassified 

3.  Title:  EVALUATION  OF  WALL  INTERFERENCE  EFFECTS  IN  A TWO- 

DIMENSIONAL  TRANSONIC  WIND  TUNNEL  BY  SUBSONIC  LINEAR 
THEORY 


4.  Personal  Author(s): 
B.  D.  Fairlie 
N.  Pollock 


5.  Document  Date: 

Written  May  1978 
Published  February  1979 


6.  Type  of  Report  and  Period  Covered: 


7.  Corporate  Authorfs):  8.  Reference  Numbers 

Aeronautical  Research  Laboratories  (a)  Task: 

(b)  Sponsoting  Agency: 

9.  Cost  Code:  54  7710 


10.  Imprint:  11.  Computer  Program(s) 

Aeronautical  Research  Laboratories,  (Title(s)  and  language(s)) : 

Melbourne 


12.  Release  Limitations  (of  the  document)  : Approved  for  public  release 


12-0.  Overseas: 


NO. 

■ 

D 

D 

■ 

B 

■ 

B 

■ 

D 

E 

13.  Announcement  Limitations  (of  the  information  on  this  page):  No  limitation 


14.  Descriptors: 

Airfoils  Linear  theory 

Transonic  flow  Wind  tunnel  tests 

Two  dimensional  flow  Wall  interference 


15.  Cosati  Codes: 
2004 
0103 


16.  ABSTRACT 

Results  of  transonic  wind  tunnel  tests  on  two-dimensional  aerofoils  are  analysed  with 
respect  to  the  effects  of  wind  tunnel  wall  interference.  The  tests  were  conducted  on  two 
geometrically  similar  models  of  each  of  two  aerofoil  sections  the  MAC  A 0012  and  the 
BGK- / sections  and  covered  a range  of  Mach  numbers  from  8-5  to  0 82  with  model 
chord  to  tunnel  height  ratios  of  0 ■ / 25  and  0 ■ 25.  Results  from  measurements  on  all  models 
in  both  solid-  and  slotted-wall  test  sections  are  corrected  for  wind  tunnel  wall  interference 
effects  by  the  application  of  classical  linearized  theory.  Tor  the  solid  wall  results,  these 
corrections  appear  to  produce  data  which  are  very  close  to  being  free  of  the  effects  of 
interference.  In  the  case  of  slotted  walls  however,  linear  theory  is  found  to  significantly 
underestimate  the  magnitude  of  blockage  interference.  The  introduction  of  an  empirical 
blockage  correction,  when  combined  with  the  linear  theory  representation  of  lift- 
interference.  enables  the  cross-flow  characteristics  of  the  slotted  walls  to  be  determined, 
and  leads  to  a successful  correction  scheme  for  the  slotted  wall  results.  Extensive  com- 
parisons with  data  from  other  sources,  both  experimental  and  theoretical,  provide  further 
verification  that  the  corrected  results  are  indeed  closely  interference  free. 


DISTR1BI  LION  I 1ST 


Copy  No. 

AUSTRALIA 
Department  of  Defence 


Central  Office 

Chief  Defence  Scientist  1 

Deputy  Chief  Defence  Scientist  2 

Superintendent,  Science  and  Technology  Programs  3 

Defence  Library  4 

Joint  Intelligence  Organization  5 

Assistant  Secretary.  DISB  6-22 

Australian  Defence  Science  and  Technology  Representative  (UK)  23 

Counsellor,  Defence  Science  (USA)  24 

Aeronautical  Research  Laboratories 

Chief  Superintendent  25 

Superintendent,  Aerodynamics  Division  26 

Divisional  File,  Aerodynamics  Division  27 

Authors:  N.  Pollock  28 

B.  D.  Fairlie  29 

Transonic  Wind  Tunnel  Group  30-35 

Library  36 

Materials  Research  Laboratories 

Library  37 

Defence  Research  Centre,  Salisbury 

Library  38 

Engineering  Development  Kstablishment 

Library  39 

RAN  Research  Laboratory 

Library  40 

Department  of  Productivity 
Government  Aircraft  Factories 

Library  41 

Statutory,  State  Authorities  and  Industry 

CSIRO  Mechanical  Engineering  Division  (Chief)  42 

SEC  Herman  Research  Laboratory  (Librarian)  Victoria  43 

Commonwealth  Aircraft  Corporation  (Manager  of  Engineering)  44 

Hawker  de  Havilland  Pty  Ltd  (Librarian),  Bankstown  45 

Hawker  de  Havilland  Pty  Ltd  (Manager),  Lidcome  46 

Universities  and  Colleges 

Adelaide  Barr  Smith  Library  47 

Professor  of  Mechanical  Engineering  48 

Australian  National  Library  49 

Flinders  Library  50 


M 


i 


d 


i 


1 


i 

I 


James  Cook  Library  51 

La  Trobe  Library  52 

Melbourne  Engineering  Library  53 

Professor  of  Mechanical  Engineering  54 

Professor  P.  N.  Joubert.  Mechanical  Engineering  55 

Monash  Library  56 

Newcastle  Library  57 

New  England  Library  58 

New  South  Wales  Physical  Sciences  Library  59 

Professor  R.  A.  A.  Bryant.  Mechanical  and  Industrial  60 

Engineering 

Queensland  Library  61 

Sydney  Professor  G.  A.  Bird,  Aeronautical  Engineering  62 

Professor  J.  W.  Roderick,  Mechanical  Engineering  63 

Professor  R.  I.  Tanner.  Mechanical  Engineering  64 

Tasmania  Engineering  Library  65 

Professor  A.  R.  Oliver,  Civil  and  Mechanical  Engineering  66 

Western  Australia  Library  67 

Professor  Alien-Williams.  Mechanical  Engineering  68 

RMIT  Library  69 

Mr.  H.  Millicer,  Aeronautical  Engineering  70 

Mr.  Pugh,  Mechanical  Engineering  71 

CANADA 

NRC  National  Aeronautics  Establishment  Library  72 

Universities 

McGill  Library  73 

Toronto  Institute  of  Aerophysics  74 

FRANCE 

AGARD  Library  75 

ONERA  Library  76 

Service  de  Documentation.  Technique  de  1’Aeronautique 

GERMANY 

ZLDI  78 

INDIA 

Defence  Ministry,  Aero  Development  Establishment,  Library  79 

Hindustan  Aeronautics  Ltd.,  Library  80 

Indian  Institute  of  Science,  Library  81 

Indian  Institute  of  Technology,  Library  82 

National  Aeronautical  Laboratory  (Director)  83 

ISRAEL 

Technion — Israel  Institute  of  Technology  (Professor  J.  Singer)  84 

ITALY 

Associazione  Italiana  di  Aeronautica  e Astronautica  (Professor  A.  Evla)  85 


i 


4 


» 4 


Ml 


JAPAN 

National  Aerospace  Laboratory,  Library  86 


Universities 

Tohuku  (Sendai)  Library  87 

Tokyo  Institute  of  Space  and  Aerospace  88 


NETIU  RI.ANDS 

Central  Organization  for  Applied  Science  Research  in  the  Netherlands  TNO,  Library  89 


National  Aerospace  Laboratory  (NLR),  Library  90 

NEW  ZEALAND 

Air  Department,  RNZAF  Aero  Documents  Section  91 

Universities 

Canterbury  Library  92 

Professor  l).  Stevenson,  Mechanical  Engineering  93 

Mr.  F.  Fahy,  Mechanical  Engineering  94 

Mr  J Stott,  Chemical  Engineering  95 

SW  EDEN 

Aeronautical  Research  Institute  96 

Chalmers  Institute  of  Technology,  Library  97 

Rungl.  Tekniska  Hogscholens  98 

SAAB.  Library  99 

SWITZERLAND 

Institute  of  Aerodynamics  (Professor  J.  Ackcret)  100 


UNITED  KINGDOM 

Aeronautical  Research  Council,  N PL  (Secretary) 

CAARC  NPl.  (Secretary) 

Royal  Aircraft  Establishment  Library.  Farnborough 

Dr  R C Lock.  RAF,  Farnborough 

Royal  Aircraft  Establishment  Library.  Bedford 

National  Physical  Laboratories  Aero  Division  (Superintendent) 

British  Library,  Science  Reference  Library 
British  Library,  Lending  Division 

Engineering  Sciences  Data  Unit  Ltd  (Mr  R H Blockley) 

Aircraft  Research  Association  l ihrary 

Science  Museum  Library 

Hawker  Siddeley  Aviation  Ltd  , Brough 

Hawker  Siddeley  Aviation  Ltd  . Greengatc 

Hawker  Siddeley  Aviation  Ltd.  kingston-upon-Thames 

Hawker  Siddeley  Dynamics  Ltd.,  Hatfield 

British  Aircraft  Corporation  (Holdings)  Ltd  . Commercial  Aviation  Division 
British  Aircraft  Corporation  (Holdings)  Ltd.,  Military  Aircraft  Division 
British  Hovercraft  Corporation  Ltd.  (East  Cowes) 

Short  Brothers 
Westland  Helicopters  Ltd. 


Universities  and  Colleges 


Bristol 


Cambridge 
Liverpool 
London 
Nottingham 
Southampton 
Strathclyde 
Cranfield  Institute  of 
Technology 
Imperial  College 


Library.  Engineering  Department 

Professor  L.  Howarth.  Engineering  Department 

library.  Engineering  Department 

Professor  J H.  Preston,  Fluid  Mechanics  Department 

Professor  A.  D Young.  Queen's  College 


Library 

Library 

Library 

Librarv 


Professor  of  Mechanical  Engineering 

Mr.  P Bradshaw.  Department  of  Aeronautics 


101 

102 

103 

104 

105 

106 

107 

108 

109 

110 
112 

113 

114 

115 

116 

117 

118 

119 

120 
121 


122 

123 

124 

125 

126 

127 

128 

129 

130 


131 

132 


UNITED  STATES  OK  AMERICA 


NASA  Scientific  and  Technical  Information  Facility 

133 

Sandia  Group  (Research  Organization) 

134 

American  Institute  of  Aeronautics  and  Astronautics 

135 

Applied  Mechanics  Review 

136 

Chemical  Abstracts 

137 

The  John  Crerar  Library 

138 

United  Technologies  Corporation,  Fluid  Dynamics  Laboratories 

139 

Battellc  Memorial  Institute,  Library 

140 

Calspan  Corporation 

141 

Universities  and  Colleges 

California 

Dr.  M Holt,  Department  of  Aero  Sciences 

142 

Cornell 

Library 

143 

Florida 

Mark  H Clarkson,  Department  of  Aeronautical  Engineering 

144 

Harvard 

Professor  A.  F.  Carrier,  Division  of  Engineering  and 

145 

Applied  Mathematics 

Johns  Hopkins 

Professor  S.  Corrsin,  Department  of  Mechanical  Engineering 

146 

Stanford 

Library,  Department  of  Aeronautics 

147 

\V  isconsin 

Memorial  Library,  Serials  Department 

148 

Brooklyn  Institute  of 

Library.  Aeronautical  Laboratories 

149 

Poly  technology 

California  Institute  of 

Library,  Guggenheim  Aeronautical  Laboratories 

150 

Technology 

Massachusetts  Institute 

Library 

151 

of  Technology 

ares 

152 

161 

